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^  A  need  exists  to  improve  the  radar  weather  display  available 
to  air  traffic  controllers.  In  support  of  this  objective,  theoretical 
and  corresponding  empirical  investigations  were  performed  to  validate 
the  feasibility  of  using  FAA  radars  (both  terminal  and  enroute)  as  the 
sensors  to  provide  the  improved  input  for  the  weather  displays.  The 
theoretical  investigation  consisted  of  a  computer  simulation  in  which 
empirical  weather  data  and  the  parameters  of  the  FAA  and  NWS  radars 
were  used  as  inputs.  By  varying  the  position  in  range  and  height  of 
reflectivity  profiles,  the  amount  of  reflectivity  that  should  be 
received  by  the  various  radars  was  computed  and  used  to  determine 
valu€!s  of  reflectivity  for  comparison  purposes.  The  empirical  program 
consisted  of  measurements  by  FAA  and  NWS  radars  of  the  power  received 
from  the  same  storms.  This  information  was  also  converted  into  reflec¬ 
tivity  values  for  comparison  purposes.  XBoth  the  theoretical  and 
empirical  results  indicate  the  FAA  rada^  (when  operated  in  a  mode 
optimized  for  weather  detection)  should  bb  capable  of  providing 
improved  input  for  use  in  the  display  of  severe  weather  in  support  of 
air  traffic  control.  NOTE;  Varying  amounts  of  equipment  modifications 
and/cr  additions  would  be  required  (depending  on  the  model  FAA  radar) 
for  simultaneously  obtaining  optimized  displays  of  aircraft  and  cali¬ 
brated  weather. 
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DETECTION  OF  SEVERE  WEATHER  BY  FAA  RADARS 


Final  Report 

EXECUTIVE  SUMMARY 

BACKGROUND 

Federal  Aviation  Administration  (FAA)  air  traffic 
controllers  need  up-to-the-minute  weather  information  to  vector 
aircraft  when  severe  storm  conditions  are  in  the  neighborhood. 
Storm  positions,  although  available  elsewhere  from  Weather 
Service  radars,  are  frequently  unavailable  in  a  timely  manner  to 
the  FAA  under  present  operating  conditions.  The  FAA  terminal  and 
enroute  radars  are  intentionally  operated  in  such  a  manner  to 
maximize  the  detection  and  tracking  of  aircraft  and  to  minimize 
the  "clutter"  of  storms.  The  FAA  recognizing  the  need  to  improve 
the  radar  weather  display  available  to  its  air  traffic  controllers 
initiated  studies  designed  to  investigate  the  feasibility  of  an 
improved  weather  display  capability  using  modified  FAA  air  traf¬ 
fic  control  (ATC)  radars. 

APPROACH  TO  PROBLEM 

Experiment 

Experiments  were  performed  by  the  FAA  to  investigate 
the  possibility  that  the  FAA  radars  might  be  reconfigured  or 
modified  to  yield  improved  weather  data  along  with  aircraft  data. 
Two  experiments  were  carried  out  in  the  summer  of  1977  in  Oklahoma 
and  Louisiana  in  which  ARSR-ID  and  ASR-8  radars  measured  storm 
intensities  simultaneously  with  Weather  Service  WSR-57  radars. 

The  WSR-57  radars  were  used  as  a  standard  against  which  the  FAA 
radars  were  compared. 

Simulation 

The  experiments  performed  in  Oklahoma  and  Louisiana 
tested  the  weather  detection  capabilities  of  two  FAA  radars.  In 
order  to  extend  the  knowledge  gained  in  the  test  to  other  FAA 
radars  and  to  guide  the  interpretation  of  the  experimental  data. 


a  computer  simulation  was  performed  that  used  empirical  weather 
data  and  the  parameters  of  the  FAA  radars  as  inputs.  The  simula¬ 
tion  involved  positioning  reflectivity  profiles  (i.e.,  reflectivity 
vs  height)  at  various  ranges  and  computing  the  amount  of  oower 
that  would  be  received  from  that  range  by  a  given  FAA  radar. 

Comparisons 

Empirical  -  WSR-57  vs  FAA 

WSR-57  and  FAA  radar  signal  levels  were  recorded  on 
\  ideo  tape  from  which  digitized  magnetic  tapes  of  the  storm  data 
V  ere  made.  These  digitized  signals  were  used  along  with  measured 
calibrations  to  compute  reflectivity  (Z)  values  in  plan  position 
indicator  (PPI)  format.  After  geometrical  colocation  of  the  data 
Irom  the  two  radars,  PPI  displays  that  were  contoured  using  the 
rational  Weather  Service  reflectivity  levels  were  compared.  In 
eddition,  comparisons  were  made  of  the  eunount  of  spatial  area 
n easured  by  each  radar  for  a  given  level. 

Simulation  vs  Empirical 

To  compare  the  empirical  results  to  the  simulation, 

WSR-57  measured  reflectivity  profiles  from  several  storms  were 
chosen  and  run  through  the  simulation  program  to  predict  what  the 
FAA  radar  should  have  measured  for  each  profile.  These  simulated 
values  were  then  compared  to  what  the  FAA  radar  actually  measured 

RESULTS 

Simulation 

Terminal 

The  simulation  results  show  that  the  ASR-8  low  beam  and 
ASR-5  radars,  whose  normal  operating  range  is  60  nmi  (111  km), 
will  cn  the  average  detect  reflectivity  levels  that  are  within  3 
dBZ  oi  the  WSR-57  when  the  maximum  gain  of  the  radars  is  posi- 
tione<l  with  the  -3  dB  point  on  the  horizon.  On  the  average,  the 
FAA  reidars  and  the  WSR-57  will  underestimate  the  maximum. 

Operating  at  ranges  less  than  60  nmi,  the  FAA  terminal  radars  are 
not  seriously  affected  by  lack  of  beaunfilling  storms  or  earth 


curvature.  The  closest  simulation  range  was  11  nmi  (20  km),  so 
that  all  results  are  for  ranges  11  nmi  and  greater.  No  attempt 
was  made  to  examine  weather  detection  at  clutter  ranges,  there¬ 
fore,  based  on  the  simulation,  no  statement  can  be  made  regarding 
the  adverse  effects  of  close-in  clutter  on  the  detection  of 
severe  weather. 

Enroute 

The  ARSR-2/3  and  ARSR-ID  enroute  radars  on  the  average 
measure  reflectivities  that  range  between  5  to  10  dBZ  down  from 
the  maximum  reflectivity  for  ranges  between  20  and  240  km  when 
the  maximum  gain  is  positioned  at  1®  elevation.  The  ARSR-2/3 
with  a  4“  (3  dB)  beamwidth  in  general  measure  reflectivities 
clciser  to  the  maximum  than  the  ARSR-lD  with  a  6“  (3  dB)  beam- 
wicth.  The  wider  beaunwidths  decrease  the  ability  to  measure  the 
ma> imum  because  as  the  storm  profile  moves  out  in  range  a  smaller 
peicentage  of  the  beam  is  filled  due  to  the  width  of  the  beam  at 
greater  ranges  and  earth  curvature.  These  becunfilling  effects 
also  cause  the  WSR-57  to  measure  less  than  the  maximum.  For  the 
maximum  gain  positioned  at  1®,  the  WSR-57  is  approximately  5  dBZ 
below  the  maximum  at  240  km. 

In  addition  to  beamfilling  and  earth  curvature  effects, 
the  effects  of  integration  within  the  beam  at  a  given  range  also 
contribute  to  the  measurement  of  reflectivity  values  by  both  the 
enroute  and  WSR-57  radars  that  are  less  than  the  maximum.  At 
most  ranges,  the  beam  contains  not  only  the  maximum  but  values 
less  than  the  maximum.  At  long  ranges  for  storm  profiles  with 
the  maximxun  at  the  surface  the  radar  may  not  illuminate  the 
maximum  because  it  is  below  the  horizon.  For  this  reason,  the 
enroute  radars  measure  values  closer  to  the  maximum  for  storms 
that  have  the  maximum  reflectivity  at  an  altitude  considerably 
abeve  the  surface. 

Positioning  of  the  enroute  radars  with  the  -3  dB  gain 
on  the  horizon  reduces  the  reflectivity  measurement,  compared  to 
1®  oositioning,  by  approximately  1  dBZ  for  the  ARSR-2  at  240  km 
and  3  dBZ  for  the  ARSR-lD  at  the  same  range.  At  close  in  ranges 


(i.e.,  less  than  80  km)  the  ARSR-ID  can  be  expected  to  measure 
reflectivities  closer  to  the  maximum  with  the  -3  dB  point  on  the 
horizon  than  it  does  with  maximum  gain  at  1“.  This  is  because 
the  storm  maximum  can  be  expected  to  be  illuminated  more  often 
with  the  higher  gain  position.  At  1®  elevation  the  maximum 
antenna  gain  is  below  the  maximum  unless  the  storm  maximum  is  on 
the  surface. 


Experimental 

Terminal 

Analysis  of  the  terminal  data  were  not  included  in  this 
report  because  of  difficulties  encountered  primarily  in  calibra¬ 
tion,  and  data  recording.  Nonetheless,  the  successful  check  of 
the  ARSR-ID  against  the  WSR-57  leads  to  the  conclusion  that 
results  of  the  simulation  of  the  terminal  radar  are  indeed  valid, 
i.e.,  that  the  ASR  radars  will  also  give  an  accurate  measure  of 
storm  intensity  beyond  the  ground  clutter  regions. 

Enroute 

Six  storm  cells  detected  by  the  ARSR-lD  radar  and  the 
WSR-57  were  analyzed  from  recorded  video  signals.  These  storms 
were  55  (30)  to  185  km  (100  nmi)  in  range  from  the  radar.  The 
storm  contours  as  seen  by  both  radars  were  found  to  be  similar  in 
shape,  location  of  the  maximum,  and  in  the  area  of  the  contoured 
levels.  On  the  average,  however,  a  difference  of  approximately  4 
dB  was  found  between  the  measured  reflectivity  values  for  the  two 
radars.  The  source  of  this  difference  has  not  been  found  but  it 
may  well  lie  in  the  calibration  of  one  or  both  of  the  radars  or 
possib.ly  in  the  positioning  of  the  maximum  gain  of  the  ARSR-lD. 
Since  this  difference  is  not  believed  to  be  fundeunental  it  can 
easily  be  taken  into  account  in  the  future  by  calibration. 

Overall  comparisons  for  the  six  storm  cells  analyzed 
from  t)ie  Oklahoma  experiment  showed  that  the  ARSR-lD  always 
measured  contour  levels  similar  to  the  WSR-57.  The  ARSR-ID  beam- 
width  appears  to  have  little  effect  on  the  measurement  of  storms 
within  185  km.  However,  the  simulation  indicated  that  measurements 


made  by  the  enroute  radars  beyond  approximately  240  km  (130  nmi) 
should  have  a  statistical  correction  applied  to  partially  correct 
for  the  effects  of  beamfilling,  earth  curvature,  and  integration. 
Such  a  correction  cannot  be  exact  because  storm  height  is  not 
known  in  real  time,  however,  it  would  enable  the  position  and 
potential  intensity  of  storms  at  long  ranges  to  be  more  accurately 
displayed  for  controller  use.  That  is,  if  a  storm  is  detected  at 
300  km  (162  nmi) ,  for  example,  a  knowledge  of  the  beamwidth  at 
that  range  would  indicate  that  the  maximum  reflectivity  in  the 
storm  must  be  above  a  certain  threshold. 

MTI 

Moving  Target  Indicator  (MTI)  data  were  analyzed  from 
the  Oklahoma  experiment.  The  analysis  indicated  that  MTI  as  it 
is  presently  configured  on  the  ARSR-ID  will  not  provide  accurate 
storm  reflectivity  estimates.  This  is  attributed  to  limiting  in 
the  MTI  receiver  because  of  the  small  dynamic  range  and  due  to 
the  unpredictability  of  storm  paths  with  respect  to  the  radar 
(i.e.,  radial  vs  tangential). 

Simulation  vs  Experiment 

Simulated  ARSR-lD  reflectivities  when  compared  to 
measured  ARSR-ID  reflectivities  gave  a  mean  difference  of  1.2  dBZ 
with  a  4.5  dBZ  deviation  about  that  mean.  A  linear  least  squared 
fit  gave  a  slope  of  .9  and  an  intercept  of  5.1.  The  general 
agreemen“  is  considered  good. 

OVERALL  f.’ONCLUSIONS 

The  following  conclusions  are  made  based  on  the  results 
from  thi:;  study. 

(1)  All  ARSR  radars  can  be  used  to  measure  and  display 

meaningful  weather  data  when  operated  in  a  mode  optimized  for 
weather  detection  within  240  km.  This  conclusion  is  based  on 
both  exporimental  results  from  the  ARSR-ID  and  simulated  results. 
The  240  )an  range  restriction  is  imposed  because  beyond  this  range 
the  effec!ts  of  earth  curvature,  partially  filled  beams,  and 
integrat  on  greatly  degrade  the  estimates  of  storm  intensity. 


(2)  All  ARSR  radars  and  the  WSR-57  will,  in  general, 
underestimate  the  maximum  reflectivity  with  the  amount  of  under¬ 
estimation  dependent  on  storm  range  and  position  of  the  maximum 
in  altitude.  For  storms  with  maximum  reflectivity  on  the  ground 
the  ur ierestimation  increases  with  increasing  3  dB  beamwidth. 

Thus,  the  WSR-57  with  the  2  degree  beam  should  be  closer  to  the 
maximum.  When  the  maximum  is  positioned  at  high  altitudes  the 
ARSR  radars  should  detect  the  storm  intensity  more  accurately 
than  the  WSR-57  radar  in  a  normal  scan  mode  (i.e.,  .5®  or  1° 
elevation  angle) .  This  result  is  significant  since  there  is 
evidence  that  severe  storms  have  maximur  intensity  at  higher 
altitudes  during  the  period  of  peak  development. 

(3)  The  ARSR  and  the  WSR-57  radars  measure  a  reflectivity 
closer  to  an  average  of  the  total  profile  reflectivity  (i.e.,  <Z>) 
than  to  the  maximum. 

(4)  The  ARSR  MTI  systems  as  they  are  presently  configured 
are  noz  suitable  for  accurate  weather  measurement. 

(5)  At  ranges  greater  than  240  km  the  ARSR  radars  should 
have  a  statistical  correction  factor  applied  to  the  received 
signal  to  account  for  partial  beamfilling,  earth  curvature  and 
integration  effects.  It  should  be  noted  that  the  WSR-57  should 
also  have  such  correction  if  operated  beyond  this  range.  The 
exact  nature  of  this  statistical  correction  would  have  to  be 
determined  based  on  a  study  of  average  heights  and  profiles  of 
severe  storms.  Since  such  a  correction  factor  would  be  statistical 
in  nature,  it  would  not  eliminate  all  reflectivity  measurement 
errors  at  long  ranges  but  rather  reduce  them. 

(6)  The  simulation  indicates  that  the  ASR-8  low  beam  and 

ASR-5  are  capable  of  detecting  weather  accurately  within  their 

normal  operating  range  and  outside  the  limits  of  ground  clutter 

(i.e.,  MTI  area).  This  statement  assumes  an  operating  mode  which 

is  suitable  for  detecting  weather  (i.e.,  linear  polarization  and 
_2 

R  STC  curve) . 

(7)  The  ASR-8  high  beam  is  not  suitable  for  weather  detection 
when  operated  as  the  sole  source  of  the  measurement.  However, 


vili 


when  it  is  used  in  a  range  gating  mode  with  the  low  beam  and 
restricted  to  close-in  ifanges,  it  may  measure  reflectivity  as 
accurately  as  the  low  beam.  The  reason  for  restriction  of  the 
high  beam  to  short  ranges  is  that  its  maximum  gain  is  positioned 
at  a  higher  elevation  angle  which  causes  it  to  illuminate  the 
upper  portions  of  storms  at  shorter  ranges  than  the  low  beam. 

The  high  beam  is  operationally  utilized  at  close-in  ranges  to 
minimize  the  reception  of  ground  clutter.  This  improves  MTI  and 
weather  reflectivity  over  ground-clutter  performance. 

(8)  The  location  and  strength  of  reflectivity  echoes  should 

aid  air  traffic  controllers  in  determining  areas  of  ir'^ense 
precipitation.  It  should  be  noted,  however,  that  re  irch  studies 
concerning  turbulence  associated  with  thunderstorms  _•  not 
established  a  clear  relationship  between  intense  rei  :tivity 
areas  and  turbulent  areas.  In  fact,  areas  of  signif  f'  .t  turbu¬ 
lence  have  been  observed  to  exist  in  the  clear  air  outside  the 
areas  of  reflectivity  associated  with  the  thunderstorm.  Based  on 
such  studies,  the  recommended  criteria  is  for  aircraft  to  avoid 
higher  reflectivity  (46  dB2  and  above  -  levels  4,  5,  and  6) 
storm  cores  by  10  to  15  miles.  Also,  it  can  be  said  that  high 
reflectivity  areas  (50  dBZ  and  above  -  levels  5  and  6)  will 
contain  intense  precipitation  and  often  aircraft  damaging  hail. 
Therefore,  if  reflectivity  data  is  properly  interpretated,  such 
displayed  information  could  materially  assist  the  air  traffic 
controllers . 

ITOTE;  Varying  amounts  of  equipment  modifications  and/or  additions 
would  be  required  (depending  on  the  model  FAA  radar)  for  simulta¬ 
neously  obtaining  optimized  displays  of  aircraft  and  calibrated 
weather . 
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DETECTION  OF  SEVERE  WEATHER  BY  FAA  RADARS 


Final  Report 


1.0  INTRODUCTION 

During  certain  seasons  in  various  parts  of  the  country 
se^'ere  storms  occur  which  do  great  damage  to  life  and  property. 
These  storms  quite  often  are  cells  which  contain  intense  precip¬ 
itation,  hail,  strong  updrafts  and  down  drafts,  and  turbulent 
wind  conditions.  The  tops  of  these  convective  storms  may  exceed 
18  km  (59  kf t)  and  several  such  cells  may  occur  aligned  in  squall 
lines  over  an  extended  range.  It  is  well  known  that  such  storms 
present  hazardous  conditions  for  aircraft  flying  in  their  vicinity; 
therefore,  the  warning  to  pilots  and  air  traffic  controllers  of 
the  exact  location  and  intensity  of  severe  storms  is  necessary 
for  aircraft  safety.  The  avoidance  of  storm  areas  by  wide  areal 
avcidance  is  not  economically  feasible.  A  more  viable  approach 
would  be  to  isolate  those  areas  with  most  probable  severe  weather, 
and  control  around  those  areas  only,  thus  reducing  time  and  fuel 
cost.  Radar,  which  has  been  used  for  this  purpose  by  the  National 
Weather  Service  (NWS)  for  many  years,  is  the  best  tool  available 
at  present  for  measuring  storm  intensity  remotely.  The  Federal 
Aviation  Administration  (FAA) ,  recognizing  the  need  to  improve 
the  radar  weather  display  available  to  its  air  traffic  controllers, 
initiated  studies  designed  to  investigate  the  feasibility  of  an 
improved  weather  display  capability  using  FAA  air  traffic  control 
(ATC)  radars. 

In  April  of  1977  the  Applied  Physics  Laboratory  (APL) 
was  given  the  task  of  evaluating,  through  theoretical  and  experi¬ 
mental  analysis,  the  detection  and  display  capability  of  the 
Federal  Aviation  Administration's  surveillance  radars.  These 
radars  have  the  sensitivity  to  detect  weather  but  are,  in  general, 
operated  in  a  mode  unsuitable  for  weather  detection.  That  is,  to 
maximize  aircraft  detection  air  traffic  controllers  usually 
operate  with  moving  target  indicator  (MTI)  and  a  sensitivity  time 
control  (STC)  curve  equivalent  to  an  r“^  range  correction  or  a 
special  adjusted  shape  that  depends  upon  siting  conditions. 
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During  periods  of  precipitation  they  may  also  operate  with  circular 
polarization.  All  of  these  modes  are  adverse  to  calibrated 
weather  detection.  In  addition  the  FAA  radars  have  wide-shaped 
antenna  beams  in  the  vertical  in  order  to  obtain  altitude 
coverage.  Wide-shaped  beams  are  not  generally  used  for  weather 
detection.  Considering  these  aspects  of  FAA  radar  operation  and 
characteristics  along  with  other  problems,  it  was  necessary  to 
analyze  their  ability  to  detect  intense  weather. 

The  problem  was  approached  in  two  parts.  The  first 
part  was  a  simulation  of  detection  using  radar  parameters  and 
empirical  weather  data  measured  by  weather  research  radars .  The 
second  part  consisted  of  two  experiments  performed  in  Oklahoma 
and  Louisiana,  respectively.  The  Oklahoma  test  was  run  to  test 
the  reflectivity  measurement  capabilities  of  an  enroute  radar  and 
the  Louisiana  test  was  run  to  test  the  terminal  radar  capabilities. 
In  both  cases  an  FAA  radar  measured  and  recorded  weather  data 
simultaneously  with  a  NWS  radar.  The  purpose  of  the  experiments 
was  to  compare  the  optimized  FAA  radar  measurement  of  weather 
reflectivity  with  that  of  the  NWS  radars  which  routinely  use 
radar  for  weather  tracking  and  forecasting. 

This  report  presents  the  results  of  both  parts  of  this 
study.  It  is  intended  that  the  results  be  used  to  deteinnine  the 
necessary  operational  and/or  physical  changes  that  must  be  made 
to  FAA  radars  for  maximizing  weather  detection.  No  evaluation  of 
these  changes  or  design  of  display  equipment  is  included  in  this 
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2.0  SIMULATION 

The  experiments  were  performed  to  test  two  types  of 
radars,  the  ARSR-ID  and  the  ASR-8.  The  tests  were  run  in  speci¬ 
fic  locations  and  thus  the  storms  were  not  necessarily  represen¬ 
tative  of  a  wide  variety  of  storms  throughout  the  country. 
Therefore,  a  simulation  study  was  designed  which  would  predict 
the  detection  capability  of  weather  by  many  FAA  radars  for  various 
empirical  reflectivity  profiles  measured  throughout  the  country. 

These  weather  detection  simulations  for  both  the  enroute 
and  terminal  radars,  using  measured  radar  data  and  real  gain 
patterns  as  inputs,  provided  insight  which  assisted  in  the  inter¬ 
pretation  of  the  results  of  the  two  experiments.  This  in  turn 
allowed  certain  deductions  to  be  made  for  other  FAA  radar  types 
operating  in  different  geographical  locations. 

2 . 1  Method  of  Simulation 

The  general  approach  to  this  simulation  was  to  use  the 
radar  equation,  radar  parameters,  and  empirical  radar  reflectivity 
profiles  to  predict  the  received  signal  from  a  storm  at  various 
radar  ranges .  Because  of  the  shape  of  the  vertical  gain  pattern 
of  the  surveillance  radars,  the  computational  procedure  will  be 
discussed  in  detail. 

Consider  the  sketch  in  Figure  2-1.  Shown  here  are  a 
series  of  elevation  angles  with  an  indicated  reflectivity  profile 
at  some  range,  R^.  Assume  that  the  profile  has  some  range  depth, 
dR,  and  values  of  as  a  function  of  altitude.  The  sketch  shows 
a  typical  shaped  beam  antenna  illumination  of  a  storm  by  indicating 
in  parenthesis  at  each  elevation  angle  the  one  way  antenna  gain 
relative  to  the  maximum  antenna  gain.  The  illumination  of  any 
elemental  volume,  v^ ,  has  associated  with  it  a  gain  depending  on 
the  height  or  elevation  angle.  This  changing  gain,  as  given  by 
the  true  antenna  pattern,  must  be  considered  in  the  simulation. 

We  designed  the  simulation  by  asking  two  questions:  The  first 
is,  what  power,  will  be  received  by  the  radar  given  that  it 
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has  detected  a  certain  profile  at  a  given  range?  Secondly,  given 
a  predicted  power  received,  what  is  the  radar  reflectivity  to  be 
associated  with  it?  We  will  proceed  with  the  equations  and 
irethods  of  computation  to  answer  these  questions. 

The  weather  radar  equation  can  be  expressed  in  several 
different  ways,  we  have  chosen  the  form  given  by  Battan  (1973) . 
Consider  the  power  received  from  an  elemental  volume. 


P  ((|>.)  = 
r  ^1 


S<*i>  >' 


where 


-  peak  transmitted  power 

R  -  slant  range  to  the  midpoint  of  the 

elemental  volume 

G,j, (4>^)  -  transmit  and  receive  gains  over  the 

elemental  volume  centered  at  angle  (J)^ 
as  given  by  the  real  measured  antenna 
pattern 

-  represents  all  system  losses  including 

■  receiver,  antenna,  and  atmospheric 

-  sum  of  all  scatterers  within  the  pulse 


volume  Vj^. 


To  simulate  what  the  FAA  radar  antenna  will  measure,  we 
will  divide  the  precipitation  volume  into  small  volumes,  v^, 
defined  by  the  horizontal  beamwidth,  0,  the  vertical  angular 
element,  A^,  and  pulsewidth,  ct/2. 

Since  we  are  calculating  the  small  volumes  as  indicated 
bi^  the  stippled  areas  in  Figure  2-1,  the  volume  is  given  by 


\  =  T  ^4-  e(f^) 
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wh  ere 


6  -  azimuthal  beamwidth 

A(t)  -  is  a  small  elevation  angle  increment 

ct/2  -  range  resolution 

T  -  pulsewidth 

In  radar  meteorology  these  scatterers  (o^)  are  in 
general  assumed  to  be  spheres  of  various  diameters  and  the  sum 
of  their  cross  sections  is  denoted  by  Z. 


L  “i  ■  7T  1*1' 


where  K  is  a  function  of  index  of  refraction  of  ice  or  water. 

'  '  n 

So  that  if  we  substitute  (2.2)  and  (2.3)  for  ^  in  Eg.  (2.1) 
we  get 


P^(^i)  = 


TT^  0  CT  |k|2  L.J.)  G^(.^.)  Gj^(4>.)  Z.  AO. 


P  (O'.)  = 
r  1 


K  G^(O'i)  Gj^(O'i)  Zi  AO. 


Where  is  the  expression  in  brackets. 

We  will  use  Eq.  (2.5)  to  compute  the  power  at  constant 
ranges  for  angles  across  the  vertical  beam.  The  total  power 
measured  by  the  shaped-beam  radar  at  any  range,  R,  is  then  given 
byr 


>.<*> '  E 


K  G^(0^)  Gj^(O'i)  Z^  AO. 


where  there  are  M  increments  of  size  A0>  used  in  the  computation. 
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This  P^(R)  value  as  computed  from  Eq.  (2.6)  represents 
the  simulated  received  power  from  the  total  pulse  volume  at  range, 
R.  We  desire  to  compute  a  reflectivity,  Z,  value  using  P^(R). 

In  order  to  do  this  let  us  write  the  radar  equation  in  an  integral 
form 

P  |K|  2  L  /* 

Pj.(R)  =  — - 3 - 5 - -  J  G,j,((t>)  Gj^(4))  Z(R)dV  (2.7) 


This  equation  is  similar  to  Eq.  (2.1)  and  is  written 
in  this  manner  in  Eq.  (2.7)  to  stress  the  Z{R)dV.  To  solve 

for  Z(R)  exactly  given  a  single  value  of  P^(R),  it  is  necessary 
to  solve  the  integral  equation  or  know  the  Z  profile  from  which 
the  P^(R)  was  measured.  In  real  life  we  will  have  no  knowledge 
of  the  real  Z  profiles^  so  an  assumption  is  made  to  compute  what 
we  call  an  "effective  reflectivity".  The  assumption  is  that  Z 
is  constant  with  height.  If  Z  is  a  constant  we  can  take  Z  out  of 
the  integral  in  (2.7)  to  obtain 


P  TT^  |k|2  L 

P  (R)  =  ^eff 

^  64  R^ 


G^(0))  Gj^((^)  dV 


(2.8) 


Note  that  should  not  be  confused  with  the  "equivalent 

reflectivity",  Z^,  usually  used  in  radar  meteorology. 

The  y  G,j,  Gj^  dV  can  be  expressed  as  a  summation 
assuming  that  gain  pattern  is  available  as  a  function  of 
elevation  angle 


f  0., 


dV  = 


(2.9) 
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So  that  (2.8)  becomes 


Pj.(R) 


'  (  nr?  ( ? 


Ac|) 


The  summation  of  gains  in  the  equation  above  are  per¬ 
formed  over  the  antenna  pattern  from  the  horizon  to  the  top  of 
the  pattern  in  elevation.  If  the  pattern  is  divided  into  N  equal 
angle  increments  the  summation  goes  from  1  to  N.  Solving  for 

ef  f 

P^(R) 

Z - - - £ -  (2.1 

K  Eg  (<^,)  G  (.^  )  H 


1 


We  have  chosen  to  use  the  integrated  gain  of  the  real 
antenna  pattern  in  the  computation  of  for  the  FAA  radars. 

Had  we  chosen  to  use  the  conventional  equation,  would  be 

computed  by 


R^  Pp(R) 


2  £n  2 


'eff 


(2.1 


Kp^MAX 


where  is  the  3  dB  elevation  beamwidth  and  2  iln  2  is  the 
Probert-Jones  correction  for  the  FAA  radar.  It  should  be  empha¬ 
sized  that  while  (2.11)  could  have  been  used  for  computation  of 
from  the  simulated  power  received,  the  integrated  gain  had 
to  be  used  in  the  computation  of  P^(R)  to  adequately  describe  a 
true  ]>ower  received.  Equations  (2.11)  and  (2.10)  will  essentially 
yield  the  same  when  the  E^.j.C'I’j^)  Gj^((()^)  A(J»  in  Eq.  (2.10)  is 

integi  ated  from  the  -3  dB  point  on  the  horizon.  This  is  not  the 
case,  however,  for  such  an  integration  when  the  -3  dB  is  not  on 
the  h<'rizon. 
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2 . 2  Parameters  for  Comparison 


Up  to  now,  only  a  simulated  has  been  discussed. 

Some  method  is  needed  to  interpret  this  Z^^^  in  terms  of  true 
storm  intensity.  To  do  this  there  must  be  something  with  which 
to  compare  the  Z^^^.  Three  quantities  have  been  chosen  for 
comparison  with  Z^^^.  For  clarity  we  will  refer  to  these  as 


(1)  Z 


WSR 


,  (2)  and  (3)  <Z'> . 


2.2.1  WSR  Reflectivity 

^WSR  obtained  by  simulating  what  a  2®  narrow  beam 
radar  would  measure  if  it  were  positioned  at  1®  elevation  angle. 
This  quantity  is  obtained  in  exactly  the  same  manner  as  the  Z^^^ 
except  that  a  2®  Gaussian  beam  pattern  and  WSR- 57  radar  parameters 
were  used. 


2.2.2  Maximum  Reflectivity 

^MAX  maximum  reflectivity  in  a  profile  and  will 

be  used  to  determine  how  well  it  is  represented  by  Z^^^. 

2.2.3  Total  Reflectivity  (<Z>) 


The  several  measures  of  Z  used  for  comparison  with  Z^^^ 


in  the  simulation  and  are  quantities  that,  in  principle 


are  observable  by  a  suitable  radar  (Weather  Service  radar  and 
maximum  reflectivity  by  a  narrow-beam  radar) .  There  is  a  heuristic 
value  in  introducing  another  Z,  which  is  independent  of  a  radar 
and  dependent  only  on  the  Z  profile.  We  call  this  "total  Z" 
(denoted  as  <2>)  ,  where  <7.>  refers  to  the  storm  profile  itself. 

This  quantity  is  defined  ns 


<Z> 


z:h)  dh 


where  h  is  height. 
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I  The  three  cuantititis  discussed  above  will  be  the 

primary  quantities  used  for  comparison  of  However,  others 

will  be  discussed  in  later  sections  as  they  arise. 

I  2.3  Computer  Program 

i 

Computer  programs  were  written  for  the  Hewlett-Packard 
9825  computer  to  perform  the  simulation  and  mean  reflectivity 
analysis.  The  main  features  of  the  programs  which  should  be 
emphasized  are:  (1)  The  simulated  power  received  is  computed  at 
.10“  elevation  angle  increments  along  the  2(h)  profile  using  the 
gain  appropriate  to  that  height,  h,  and  summed  to  the  top  of  the 
storm,  (2)  The  simulated  power  received  from  the  FAA  radar  is 
turned  into  a  using  Eq.  (2.10),  (3)  Z^^^^  is  computed  in  a 

manner  similar  to  Z^^^  using  the  WSR-57  gain  pattern,  (4)  Each 
profile  is  positioned  in  range  at  20  km  increments. 

The  simulation  program  computed  the  three  variables, 
^MAX'  ^WSR'  ^eff'  profiles  and  computed  the  difference 

between  the  first  two  and  Another  program  was  then  used 

to  compute  the  mean  difference  between  these  parameters  and  their 
correlation. 

2 . 4  The  FAA  Radars 

Six  radar  types  were  considered  in  the  simulation.  The 
nominal  radar  parameters  for  these  radars  are  shown  in  Table  1. 
The  antenna  gain  patterns  for  each  of  the  radars  is  given  in 
Appendix  B,  The  main  characteristics  of  all  FAA  radars  are 
their  wide-shaped  ver-c-cai  beams  and  narrow  Gaussian  horizontal 
beams.  They  all  utilize  MTI  for  ground  clutter  cancellation  and 
may  employ  various  othrr  weather  clutter  rejection  methods  such 
as  circular  polarizanacr  C7AR,  etc.  In  the  case  of  the  ASR-8 
and  ARSR-3,  there  are  r'  o  beams  on  the  radars  called  high  and  low 
beams.  During  .ormai  o,  -c.  cions  the  high  beam  and  low  beam  are 
often  used  in  a  I'ange  mode.  Transmission  occurs  via  the 

low  .intepna  bean,  anc  .-Ccpcicn  is  achieved  on  either  the  low  or 
high  beam. 


ARSR-2/1E  ARSR-3  ASR-8  ASR-4/5/6/7  WSR-57 
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The  gain  patterns  for  each  of  these  radars  were  read 
at  .5°  elevation  angle  increments  and  recorded  on  HP  9825  cassette 
tapes  for  use  in  the  program. 

2 . 5  Empirical  Weather  Data 

Many  research  groups  have  been  involved  in  the  study  of 
thunderstorms  using  radar  as  a  remote  sensing  tool  to  measure  the 
intensity  and  spatial  distribution  of  precipitation  within  storms. 
For  this  simulation,  we  have  chosen  profiles  of  radar  reflectivity 
vs  altitude  from  many  of  these  researchers. 

There  were  37  profiles  used  in  the  simulation  and  the 
set  contains  basically  three  types  such  as  those  shown  in  Figure 
2~2.  They  are:  (A)  The  reflectivity  is  constant  from  the  ground 
to  som.e  altitude  and  decreases  with  altitude  from  that  point  to 
the  top  of  the  storm.  (B)  The  reflectivity  is  lower  on  the 
grounl  than  it  is  aloft.  (C)  The  reflectivity  is  maximum  on  the 
ground  and  decreases  gradually  to  the  top  of  the  storm.  The 
shape  of  these  profiles  may  be  indicative  of  the  type  of  storm 
although  little  research  has  been  done  to  determine  this. 

Some  of  the  profiles  are  of  thunderstorms  with  hail. 

The  size  that  hail  reaches  can  be  much  greater  than  raindrop 
sizes  .  Such  hail  with  a  coating  of  water  will  produce  an 
especially  high  return  (reflectivity  is  proportional  to  the  sixth 
power  of  a  linear  dimension) .  The  profiles  used  in  the  simula¬ 
tion  were  taken  from  thunderstorms  for  che  most  part  because  it 
is  pi imarily  these  storms  that  are  potentially  dangerous  to 
aircraft. 

2 . 6  Discussion  of  Simulation  Results 

In  the  simulation,  37  reflectivity  profiles  were 
anali  zf  .  as  mentioned  p'reviously.  In  this  section  the  three 
showr  in  Figure  2-2  wi]l  be  used  for  discussion  along  with  the 
ARSR- 2  and  ASR-8  radars.  Profile  A  was  positioned  in  range  and 
a  sinulated  detection  by  the  ARSR-2  and  ASR-8  radars  was  computed. 
In  acdition,  the  paramt'ters  and  <Z>  were  computed. 
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Reflectivity,  as  a  function  of  range,  is  shown  for 
profiles  A,  B,  and  C  and  the  various  radars  in  Figures  2-3  through 
2-8.  Two  radar  positions  are  shown  for  the  ARSR-2.  One  is  with 
the  position  of  maximxim  gain  (PMG)  at  1“  and  the  other  with  PMG 
such  that  the  -3  dB  point  of  the  gain  pattern  is  on  the  horizon, 
PMG  at  2®.  The  ARSR-1  is  shown  for  a  1®  PMG. 

2.6.1  Individual  Profiles  -  Terminal 

The  values  of  the  reflectivity  calculated  for  the  ASR-8 
low  becun  at  1®  PMG,  using  profile  A,  lie  within  3  dBZ  of  the  max- 
imur  at  all  ranges  and  are  within  2  dBZ  of  <Z>,  as  can  be  seen  in 
Figure  2-3.  The  ASR-8  low  beam  positioned  at  2.5®  PMG  measured  a 
reflectivity  within  4  dBZ  of  the  maximum  at  all  ranges.  The 
WSR-57  measures  the  profile  within  1  dBZ  at  all  ranges. 

The  3  dB  beamwidth  of  the  radars  and  the  portion  of  the 
vertical  profile  being  illuminated  determines  the  radar's  estimate 
of  reflectivity.  Profile  A  is  constant  from  the  ground  to  4.6  km 
(15  kft)  and  decreases  at  3  dBZ  per  km  from  this  height  to  13  km 
(43  kft).  Therefore,  the  WSR-57  (1®  PMG)  is  illuminating  only 
the  maximum  at  ranges  less  than  111  km  (60  nmi) .  The  ASR-8  low 
beam  (1®  PMG)  starts  to  illuminate  the  upper  portion  of  the  storm 
at  approximately  69  km  (37  nmi)  and  therefore  measures  a  reflect¬ 
ivity  lower  than  the  maximum  beyond  this  range.  As  the  storm 
moves  out  in  range  more  of  the  upper  portion  is  illuminated  and 
the  ASR-8  low  beam  reflectivity  decreases  further  from  the  maximum 
The  ASR-8  high  beam  (4.5®  PMG)  starts  to  drop  below  the  maximum 
at  an  even  shorter  range  for  the  same  reason.  In  addition,  the 
ASR-8  high  beam  (4.5®  PMG)  hes  only  a  portion  of  the  beeun  filled 
starting  at  approximately  lOC  km  (54  nmi)  range  and  so  for  this 
radar  at  4.5®  PMG  the  reflectivity  is  also  being  reduced  due  to 
partial  beamfilLing.  The  sane  thing  applies  to  the  ASR-8  high 
beam  ;6"  PMG)  with  a  greater  effect  on  the  reflectivity  measure. 

Good  overall  -eflectivity  performance  is  achievable  by 
I  ange  gating  of  the  AS.-.  6  high  and  low  beams  as  a  function  of 


Figure  2-3. 


Simulated  reflectivity  for  ASR-8  and 
WSR-57  -  Profile  A. 


Page  Sixteen 


azimuth  (RAG)  such  that  the  high  beam  is  only  utilized  at  close-in 
ranges  when  necessary  to  minimize  the  reception  of  ground  clutter. 
This  permits  improvements  in  MTI  and  weather  reflectivity 
performance. 

The  results  from  profile  B,  shown  in  Figure  2-4,  show  a 
different  effect.  In  this  case  the  profile  has  its  maximum 
reflectivity  at  6.2  km  (20  kft)  height.  The  terminal  radars  with 
wider  beams  illuminate  the  maximum  at  shorter  ranges  than  the 
WSR-57  with  its  2®  beeim.  However,  because  the  reflectivity  is 
less  than  the  maximum  below  6  km  the  reflectivity  averaged  over 
the  beam  is  less  than  the  maximum  for  all  radars  at  all  ranges. 

The  ASR-8  high  beam  (6®  PMG)  falls  off  from  the  maximum  more  at 
120  km  because  the  beam  is  partially  filled. 

The  reflectivities  for  profile  C  shown  in  Figure  2-5 
have  a  trend  very  similar  to  those  of  profile  A.  This  should  be 
expected  since  profile  C  has  its  maxim\am  on  the  ground  and  the 
maximvun  is  illuminated  starting  at  very  short  ranges. 

It  should  be  noted  that  for  all  three  profiles,  <Z>  is 
less  than  the  maximum  because  it  is  an  (ffective  average  over  the 
total  profile.  In  addition,  values  of  reflectivity  for  the 
various  radars  are  closer  to  <Z>  than  the  maximum.  Notable 
exceptions  exist  due  to  the  particular  spatial  relationships  for 
profile  A  betv/een  the  storm's  maximum,  the  points  of  the  maximum 
antenna  gains,  and  the  antenna  vertical  beeun  shape.  In  general, 
radars  do  not  illuminate  profiles  that  are  constant  with  respect 
to  reflectivity  within  the  beam  (as  is  the  case  for  profile  A) , 
hence  all  radars  including  the  WSR-57  measure  values  of  reflectivity 
closer  to  <Z>  than  to  the  maximum, 

2.6.2  Individual  Prof iles  -  Enroute 

The  FAA  enroute  radars  and  the  WSR-57  are  affected  by 
part-al  beamfilling  and  earth  curvature.  The  results  in  Figures 
2-6  through  2-8  show  this  effect  rather  vividly.  Because  these 
effects  cause  increasing  error  with  increasing  range  on  all  three 
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profiles,  we  will  restrict  discussion  of  the  results  to  ranges 
less  than  approximately  240  km  (130  nmi) .  This  240  km  range  is 
close  to  the  normal  operating  range  of  the  WSR-57. 

Consider  the  results  of  profile  A  in  Figure  2-6. 

Inside  of  120  km  (64.8  nmi)  the  ARSR-2  (1“  and  2®  PMG)  is  within 
about  4  dBZ  of  the  maximum.  Beyond  this  range  the  ARSR-2  reflec¬ 
tivity  at  2®  PMG  decreases  more  rapidly  than  the  AprR-2  at  1®  PMG 
The  ARSR-ID  at  1®  PMG  has  a  steady  fall  off  from  the  maximum  out 
to  240  km  and  measures  reflectivities  comparable  to  the  ARSR-2 
(2®  PMG)  at  240  km.  The  WSR-57  also  falls  off  from  the  maximum 
starting  at  approxir.ately  140  km. 

The  relative  effects  of  illuminating  the  total  profile 
A,  and  beeunfilling  can  be  illustrated  by  considering  the  reflec¬ 
tivity  measurements  at  200  km  (107.9  nmi)  in  Figure  2-6.  The 
ARSR-ID  (1®  PMG)  is  about  75%  filled  and  is  illuminating  the  pro¬ 
file  from  about  2.3  km  (7,5  kft)  to  the  top.  The  reflectivity 
measured  by  the  ARSR-ID  (1*  PMG)  at  200  km  is  56  dBZ  or  approxi¬ 
mately  6.5  dBZ  down  from  the  maximum.  About  1.2  dBZ  of  this 
difference  can  be  attributed  to  partial  beamfilling  while  the 
rest  is  due  to  averaging  of  the  profile  within  the  beam  above 
the  horizon.  Similar  reasoning  can  be  used  for  the  behavior  of 
ARSR-2  reflectivity  measures. 

Profile  B  results  are  shown  in  Figure  2-7.  As  with  the 
terminal  radars,  the  wider  besuns  of  the  enroute  radars  illuminate 
the  maximum  at  shorter  ranges  than  the  WSR-57.  However,  as  men¬ 
tioned  in  previous  discussion  of  this  profile  (i.e.,  terminal 
case) ,  all  radars  fall  short  of  measuring  the  maximum  at  all 
ranges.  Regardless  of  the  effects  of  earth  curvature  or  partial 
boamfilling,  there  is  no  range  inside  of  240  km  where  only  the 
maximum  value  is  within  the  beam.  Therefore  an  averaged  quantity 
is  measured.  The  AP9R-2  (1®  PMG)  measures  values  7  dBZ  down 
from  the  maximum  at  40  km  as  compared  to  8  dBZ  down  by  the  WSR-57 
The  ARSR-2  (1®  PMG)  reflectivities  increase  from  this  range  up  to 


e  maximum  (Zmax) 


tnulated  reflectivity  for  ARSR-ID,  ARSR- 
id  WSR-57  -  Profile  A. 
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a  value  which  is  5  dBZ  below  the  maximum  at  160  km  and  decrease 
from  this  point  on.  The  ARSR-2  (2®  PMG)  measures  within  3.5  dBZ 
of  maximum  at  120  km.  The  ARSR-ID  (1°  PMG)  measures  within  4  dBZ 
of  maximum  at  120-160  km.  The  range  at  which  the  particular 
radar's  reflectivity  is  closest  to  the  maximum  varies  depending 
on  the  width  and  position  of  the  beam.  The  WSR-57  measures  closest 
to  the  maximum  at  200  km  where  it  is  within  2  dBZ  because  its 
full  narrower  beamwidth  illuminates  the  best  portion  of  the  profile 
at  that  range. 

The  results  of  profile  C  shown  in  Figure  2-8  indicate 
that  none  of  the  radars  measure  the  maximum  exactly  at  any  range. 
Again,  this  is  because  an  average  of  the  total  profile  within 
the  beam  is  being  measured  rather  than  the  maximum.  This  is  evi¬ 
denced  by  the  fact  that  the  radars  measure  values  closer  to  <Z> 
than  to  the  maximum.  The  general  trend  in  reflectivity  values 
from  profile  C  is  a  decrease  with  range  attributed  to  the  portion 
of  the  profile  being  illuminated,  partial  beamfilling  and  earth 
curvature.  The  amount  of  the  decrease  depends  on  the  range  and 
the  particular  radar  making  the  measurement. 

Results  from  the  ARSR-3  have  not  been  included  in 
Figures  2-6  through  2-8  because  the  values  are,  in  general, 

within  1  dBZ  of  the  ARSR-2  at  the  same  PMG. 

2.6.3  Mean  Reflectivity  for  Profile  Set 

The  data  chosen  for  the  profile  set  used  in  this  study 
refle::ts  research  data  from  the  various  parts  of  the  United 
States.  It  does  not  represent  a  statistical  sample.  Storm 
research  using  radars  is  still  in  the  early  stages  with  respect 
to  a  statistical  model  or  models  representative  of  various  areas 
of  the  country.  Therefore,  the  reason  for  discussing  the  mean 
refl»  ctivity  for  37  profiles  is  not  for  universal  application  of 
the  I esults  in  a  statistical  sense.  Rather  the  presentation  of 
the  .  aan  for  the  profile  set  allows  a  convenient  method  of  dis¬ 
cussing  average  results  a  given  radar  would  produce.  The  profile 
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set  does  include  a  wide  range  of  reflectivities,  shapes,  and 
rax ir urn  altitudes.  However,  none  of  the  profiles  have  maximum 
reflectivities  less  than  30  dBZ  and  it  is  highly  improbable  that 
a  stcrm  cell  with  maximum  reflectivity  less  than  30  dBZ  is  severe. 


The  procedure  for  obtaining  the  mean  reflectivities 
discussed  in  this  section  is  as  follows:  (1)  results  such  as 
those  shown  in  Figure  2-6  were  computed  for  each  of  the  thirty- 


seven  profiles.  That  is,  a  and  Z^^^^  were  computed  as  a 


function  of  range  for  each  profile.  (2)  A  Z  and  <Z>  were 


obtained  for  each  profile.  (3)  The  mean  reflectivities  were  then 
obtained  by  averaging  the  37  Z^^^  values  at  each  range  and  the 
values  at  each  range.  (4)  <Z>  and  are  quantities 


37  Z 


WSR 


obtained  from  averaging  the  37  values  of  <Z>  and  respectively 


To  interpret  Z^^^^  in  terms  of  the  radars  one  can  say  that  if  a 


given  radar  at  some  range  had  measured  Z  exactly  for  all  37 

MAX 


profiles,  then  its  mean  reflectivity  would  be  equal  to  Z 


MAX' 


Therefore,  comparing  the  mean  reflectivities  of  the  data  to  Z 


MAX 


gives  an  indication  of  the  average  performance  of  the  radar  in 
measuring  maximum  reflectivities  for  the  profile  set. 


The  ASR-8  Low  and  ASR-S  have  very  similar  gain  patterns 
so  it  should  be  expected  that  the  results  would  be  similar.  The 
data  show  (see  Figs.  2-9  and  2-10)  that  on  the  average  for  the 
profile  set  the  ASR-8  low  and  ASR-5  (outside  the  MTI  area)  are 
within  1  dBZ  of  <Z>  at  all  ranges  for  1°  PMG  and  within  1.5  dBZ 
of  <Z>  for  2.5°  PMG  (-3  dB  on  horizon) .  On  the  average  these  two 
radars  underestimate  Z^^  by  between  3  and  4  dBZ  at  1°  PMG  and 
up  to  5  dBZ  at  2.5°  PMG  at  terminal  ranges.  The  WSR- 57  under¬ 
estimates  .  by  between  2  and  3  dBZ  and  is  within  1.5  dBZ  of 
<Z>  for  1°  PMG,  It  should  be  noted  that  a  1°  PMG  for  the  WSR-57 
is  the  same  as  the  -3  dB  point  on  the  horizon. 


The  ASR-8  high  beam  is  positioned  such  that  at  both 
antenna  tilt  positions  analyzed  there  is  a  greater  decrease  in 
reflectivity  with  range  than  for  the  corresponding  ASR-8  low 
beam.  The  ASR-8  low  beam  provides  a  much  better  reflectivity 
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measure  than  the  ASR-8  high  beaun  at  all  ranges  greater  than  60  km 
(32.4  nmi)  .  Refer  to  Section  2.6.1  for  range  gating  of  ASR-8 
beams  to  achieve  overall  improvements  in  the  MTI  and  weather 
performance . 

The  mean  reflectivities  which  are  shown  in  Figures 
2-11  and  2-12  for  the  enroute  radars  reflect  the  influence  of 
the  3  dB  beeimwidth  as  a  function  of  range.  At  ranges  less  than 
about  80  km  (43.2  nmi)  the  ARSR-ID  with  its  6®  beamwidth  illumi¬ 
nates  a  larger  portion  of  the  storms  than  the  WSR-57  and  ARSP-2 
at  both  beam  positions  considered.  This  larger  beam  is  advan¬ 
tageous  in  those  cases  where  the  storms  peak  at  altitude.  Inside 
80  km,  this  is  reflected  in  the  ARSR-ID  reflectivity  values  which 
are  closer  to  „  than  the  other  two  radars.  For  example,  with 
the  -3  dB  point  on  the  horizon  (Figure  2-12)  at  80  km  the  differ¬ 
ences  between  „  for  the  ARSR-ID,  ARSR-2,  and  WSR-57  are 
respectively,  1.2,  3.0,  and  4.0  dBZ.  As  the  range  increases  the 
ARSR-2  and  the  WSR-57  measure  closer  to  the  than  the  ARSR-ID 

with  the  estimate  of  Z^^  degrading  increasingly  with  range  for 
all  the  radars.  At  240’km  (129.6  nmi),  the  ARSR-lD,  ARSR-2,  and 
WSR-57  are  on  the  average  12,  8,  and  4.5  dBZ  low  in  the  measure- 
ment  of 

The  data  in  Figure  2-12  show  that  on  the  average  for 
the  profile  set,  the  ARSR-ID  at  3®  PMG  is  12  dBZ  down  from  Z^^^, 
and  the  ARSR-2  at  2®  PMG  is  8  dBZ  down  from  at  24  0  km,  while 

the  VJSR-57  at  1®  PMG  is  4.5  dBZ  down  from  Z„- „  at  240  km.  The 
ARSR-2  at  2®  PMG  is  about  4  dBZ  down  and  the  ARSR-ID  at  3®  PMG 
is  about  8  dBZ  down  from  <Z >  at  240  Vm  range  as  compared  to  only 
1  dBZ  down  fror.  <”  >  at  240  km  by  the  WSR-57  at  1®  PMG.  This  may 
not  be  as  serious  as  it  sounds,  because  knowing  that  the  potential 
Cor  this  exists,  average  corrections  vs  range  can  be  made.  In 
othei  words,  unlike  r eteorological  research  where  accurate  dBZ 
measures  are  needed  lor  computation  of  rain  rate  and  the  dynamics 
of  atmospheric  physics,  the  controllers  require  a  knowledge  of 
location  and  potenticl  severity  of  storms. 


Figure  2-11.  Mean  reflectivity  vs  range  for  ARSR-ID, 
ARSR-2,  and  WSR-57. 


Figure  2-12.  Mean  reflectivity  vs  range  for  ARSR-lD,  ARSR-2,  and  WSR-57. 

Maximum  gain  positioned  with  -3  dB  point  on  the  horizon. 
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The  ARSR-3  results  were  not  included  in  Figures  2-11 
and  2-12  because  they  are  within  0.5  dBZ  of  the  ARSR-2  at  all 
ranges . 

It  should  be  noted  that  the  simulated  power  received  was 
calculated  based  on  the  direct  power  only  that  would  be  received 
for  the  selected  antenna  tilt  angle  (i.e.,  no  consideration  for 
ground  reflected  energy),  was  calculated  based  on  the  free 

space  antenna  pattern  for  the  selected  tilt  angle  (i.e.,  integra¬ 
tion  of  the  power  above  the  horizon  without  regard  to  ground 
reflections) . 

It  was  recognized  that  both  r  .  (calculated)  and 
(calculated)  for  the  simulation  would  be  altered  to  some  extent 
because  of  the  eimount  of  energy  added  to  the  pattern  above  the 
horizon  due  to  ground  reflections.  The  magnitude  of  the  reflected 
energy  at  each  elevation  angle  above  the  horizon  is  a  function  of 
antenna  height,  wavelength,  and  coefficient  of  reflection  of  the 
ground.  In  order  to  establish  bounds  on  this  effect  a  separate 
worst  case  analysis  was  performed  on  the  effect  of  ground  reflections 
assuming  a  coefficient  of  reflection  (p)  equal  to  1.  The  results 
of  this  analysis  are  discussed  in  Section  5.2  in  conjunction  with 
the  experimental  data  but  the  results  could  be  applied  to  the  sim¬ 
ulation  as  well.  The  discussion  of  the  analysis  is  delayed  until 
later.  Based  on  this  analysis,  a  general  statement  can  be  made 
that  if  the  worse  case  (i.e.,  p  =  1)  applies,  the  ARSR-lD  simu¬ 
lated  reflectivities  for  a  1"  PMG  would  be  increased  by  2  dBZ  on 
the  average  between  40  )an  (21.6  nroi)  and  160  km  (86  nmi)  .  The 
WSR-57  simulated  reflectivities  would  not  be  increased  at  normal 
operational  PMG's  (above  0.5**). 
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3.0  EXPERIMENTS 

3 . 1  New  Orleans  -  Slidell 

During  August  of  1977  an  experiment  was  conducted  in 
Louisiana  to  compare  simultaneous  weather  detection  by  the  FAA 
ASR-8  terminal  radar  at  New  Orleans  International  Airport  (NOIA) 
and  the  NWS  WSR-57  at  Slidell,  Louisiana.  The  relative  location 
of  the  two  radars  is  shown  in  Figure  3-1 (a) . 

A  video  tape  recorder  was  used  at  the  WSR-57  site  to 
record  raw  log  video  during  a  data  run.  A  data  run  consisted  of 
the  antenna  being  scanned  in  a  PPI  mode  at  each  of  the  following 
elevations  to  permit  definition  of  vertical  profiles;  1®,  5®, 

9®,  13®,  11®,  7®,  and  3®.  At  the  beginning  of  each  day  a  radar 
calibration  was  recorded  on  video  tape  at  3  dB  power  intervals 
and  after  each  run,  a  short  calibration  of  four  power  levels  was 
recorded.  No  radar  modifications  were  made  at  the  WSR-57  site 
and  the  radar  was  operated  exactly  as  the  NWS  operates  it.  The 
WSR-57  radar  was  available  for  the  experiment  only  five  minutes 
out  of  each  half  hour  because  it  was  an  NWS  operational  facility. 

At  the  ASR-8  site  certain  modifications  were  necessary 
in  the  receiver  to  detect  orthogonal  circular  polarization  on  the 
standby  channel  of  the  radar,  since  the  radar,  in  practical  oper¬ 
ation,  often  uses  circular  polarization  which  attenuates  weather 
signals.  This  modification  was  the  only  major  one;  however, 

equipment  was  built  to  switch  modes  of  the  standby  radar  channel 

-2 

automatically  during  a  data  run  and  the  R  STC  curve  was  used 
instead  of  the  usual  R  .  An  automatic  switching  and  recording 
sequence  was  necessary  for  the  standby  channel  at  NOIA  to  enable 
data  collections  which  would  later  permit  examination  of  various 
methods  of  optimizing  the  7\SR-8  system  for  weather  detection  for 
comparison  with  the  five  minute  WSR-57  data.  A  data  run  for  the 
ASR-8  consisted  of  continual  video  tape  recording  of  the  log 
normal  video  and  MTI  video  throughout  a  WSR-57  five  minute  data 
run,  where  the  beginning  of  a  run  was  signaled  vocally  via 
telephone. 
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True  north 


New  Orleans 


True  north 


Norman 


(b)  Relative  location  of  radars  in  Oklahoma  test 


Figure  3-1.  Relative  locations  of  radars  in  Oklahoma 
and  Louisiana  experiments. 
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The  ASR-8  has  two  feed  horns.  The  low  beam  is  used  for 
transmit  and  receive  and  the  high  beam  is  used  for  receive  only. 

An  automatic  sequence  consisted  of  receiving  four  360“  scans  of 
the  low  beeun  with  no  signal  attenuation;  four  scans  of  low  beeim 
with  the  signal  attenuated;  four  scans  of  high  beam  with  no  signal 
attenuation;  and  four  scans  of  high  beam  with  attenuation.  The 
purpose  for  recording  low  beeim  and  high  beam  in  sequence  was  to 
investigate  the  feasibility  of  using  the  high  and  low  beams  in  a 
range  gating  mode  to  aid  in  detection  of  weather  over  ground 
clutter.  That  is,  if  high  beeuti  detection  at  close  ranges  were 
sufficiently  accurate,  it  could  be  used  over  ground  clutter  areas 
where  its  higher  beam  position  reduces  or  avoids  clutter.  A 
switch  would  then  be  made  to  the  low  beam  at  ranges  outside 
ground  clutter.  The  attenuation  was  placed  in  the  sequence  to 
investigate  the  use  of  attenuation  in  more  accurately  measuring 
reflectivity  values  with  MTI  particularly  when  received  with 
ground  clutter.  Since  the  MTI  receiver  dynamic  range  is  small, 
higher  reflectivities  and  weather  returns  with  strong  ground 
clutter  saturate  the  MTI  receiver.  It  was  hoped  that  attenuation 
would  bring  higher  signals  out  of  saturation  and  permit  more 
accurate  measurements  of  weather  reflectivity.  On  some  occasions 
this  automatic  sequence  was  replaced  in  a  data  run  with  other 
changes  of  standby  radar  channel  operation  to  collect  data  for 
analysis  concerning  the  effect  of  specific  radar  circuitry  on 
weather  detection  and  calibration,  such  as  MTI  velocity  response. 

Calibration  levels  were  recorded  for  the  ASR-8  in  a 
manner  similar  to  the  WSR-57. 

Data  runs  were  recorded  on  several  days  during  the 
test.  At  the  end  of  the  experiment  the  data  were  carried  to  NAFEC 
to  be  digitized  and  sent  to  APL  for  processing  and  analysis. 

3 . 2  Oklahoma  City  -  Norman 
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detection  of  weather  to  the  WSR-57,  located  at  Norman.  The  rela¬ 
tive  locations  of  these  radars  is  shown  in  Figure  3-1 (b).  The 
experiment  is  described  in  detail  in  Appendix  C. 

In  an  operation  similar  to  the  Louisiana  experiment  a 
video  tape  recorder  was  used  to  record  MTI  and  raw  log  normal 
video  signal  at  the  ARSR-lD  site.  A  video  tape  recorder  was 
also  used  at  the  WSR-57  site  to  record  raw  log  video.  Antenna 
elevation  information  was  also  recorded  along  with  radar  cali¬ 
bration  signals.  A  data  run  consisted  of  scanning  the  WSR-57 
in  stepped  elevation  sequences  while  the  ARSR  scanned  continuously 
in  various  operational  configurations.  The  elevation  sequences 
at  Norman  depended  on  the  range  of  the  storm  and  were  either  0®, 
2®,  4®,  etc.,  to  the  top  of  the  storm  or  1®,  5®,  9®,. ..,7®,  3®, 

0®.  No  modifications  were  made  to  WSR-57  and  it  was  used  in  its 
usual  operational  mode. 

The  ARSR-ID  at  Oklehoma  City  is  a  radar  normally  used 
for  teaching  so  that  at  certain  times  it  could  be  dedicated  to 
the  experiment.  No  automatic  sequence  of  modes  was  run  in  this 
test.  Rather  various  radar  configurations  were  run  depending  on 
the  rain  conditions  and  the  configurations  to  be  tested.  The 
configurations  for  which  data  were  recorded  are  shown  in  Table  2. 
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NOTE:  These  configurations  were  taken  from  Table  1  of  Appendix  C.  The  configuration  numbers 
are  not  sequential  because  some  planned  tests  were  eliminated  before  the  experiment 
began. 
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4 . 0  DATA  REDUCTION  -  ENROUTE 

Video  was  recorded  for  simultaneous  observations  ol 
storms  by  the  ARSR-ID  and  the  WSR-57  in  order  to  compare  the 
reflectivities  measured  by  each  radar.  Power  levels  along  with 
azimuth  and  elevation  information  were  digitized  and  recorded  on 
tape  for  six  selected  storm  cells.  These  data  were  then  processed 
via  a  computer  program.  The  method  of  data  reduction  will  be 
discussed  here. 

4 . 1  Digitizing  Raw  Data 

The  data  on  both  the  ARSR  tapes  and  the  WSR-57  tapes 
were  digitized  at  range  intervals  of  .869  usee  and  at  the  indi¬ 
vidual  pulse  repetition  frequencies.  The  ARSR  has  a  prf  of  360 
with  a  6  rpm  scan  rate  which  results  in  a  digitized  data  resolu¬ 
tion  of  130  meters  in  range  and  .1“  in  azimuth.  The  corresponding 
resolution  for  the  WSR-57  is  130  meters  (.869  ysec)  and  .11®  since 
it  has  a  long  pulse  prf  of  164  and  scan  rate  of  3  rpm.  The  equip¬ 
ment  used  to  digitize  the  ARSR-ID  and  WSR-57  data  is  limited  in 
the  number  of  data  points  that  can  be  digitized.  Therefore,  the 
size  of  the  storm  cell  digitized  varies  depending  on  its  position. 
Those  stormt  where  the  majority  of  the  cell  could  be  digitized 
were  digitised  provided  they  met  all  other  criteria  in  terms  of 
intensity,  1  eight,  etc.  The  sizes  of  cells  digitized  can  be  seen 
in  the  figu:  es  of  Section  5. 

4.2  O tlibration  -  ARSR-1  and  WSR-57 

Calibration  curves  were  digitized  for  the  ARSR  and 
WSR-57.  Th(?  power  levels  recorded  on  other  days  were  then 
checked  fron  the  analog  signal  to  be  certain  no  significant 
changes  occurred  from  day  to  day.  These  calibration  curves  for 
mean  power  levels  averaged  over  the  entire  digitized  signal  are 
shown  plotted  in  Figure  4-1.  The  ARSR  calibration  was  recorded 
with  the  STC  curve  in  the  system,  therefore,  the  power  levels 
were  digitized  at  a  ra age  beyond  the  influence  of  the  STC  curve. 


WSR  57 
ARSR  ID 
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Figur  4-1.  Calibration  curves  for  ARSR-ID  and  WSR-57 
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The  WSR-57  was  recorded  before  injection  of  the  STC  signal,  so 
that  it  was  not  a  factor.  The  mean  noise  level  which  was  obtained 
by  averaging  a  digitized  noise  sample  from  the  same  time  period 
was  removed  from  both  the  curves.  The  received  pov^/er  from  the 
sun  was  recorded  on  video  tape  for  each  of  the  radars  as  an  addi¬ 
tional  calibration  check  but  the  signals  have  not  been  digitized 
to  date. 

4 . 3  Averaging 

The  computer  program  was  designed  to  read  a  scan  of 
ARSR  data  from  the  tape  and  average  I  samples  in  range  and  J 
samples  in  azimuth,  where  I  and  J  were  inputs  to  the  program  and 
varied.  The;  averaging  was  done  in  tape  units  (i.e.,  digitized 
values  between  0  and  1024)  and  a  statistical  correction  factor 
was  applied  to  the  radar  equation  to  correct  for  log  averaging. 
Similarly,  a  scan  from  the  WSR  was  read  from  a  tape  and  averaged 
for  I  and  J  samples.  It  was  these  averaged  values  that  were 
used  for  computing  reflectivity  and  comparing  areas  of  colocated 
reflectivity.  For  the  WSR-57  the  samples  were  averaged  in 
elevation  and  range  for  one  portion  of  the  analysis. 

4 . 4  Computation  of  Reflectivity 

The  simulation  section  of  this  report  discussed  the  form 

of  the  radar  equations  used  for  the  ARSR  in  this  analysis.  This 

equation  using  the  integrated  gain  is  given  below  (conversion 

18  6  6 

factors:  1852  meters  =  1  nautical  mile;  10  mm  =  1  m  ); 


24  2 

,  3.43  X  10  R  P 

Z  ,-(mmVm'’)  =  - ^ - -  (4.1) 


IK 


CT 

T" 


ARSR 


256  \- 


where  P^  and  P^  are  in  watts;  0  is  in  radians;  x  is  in  seconds; 
A  is  in  meters;  and  R  is  in  nautical  miles. 
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Table  3  lists  the  parameter  values  for  the  two  radars. 

The  ARSR  system  losses  are  included  in  the  calibration  and  the 
onli  loss  included  in  the  computation  is  -2.54  dB  for  log  averaging. 
Included  also  is  a  correction  of  .7  dB  for  the  Gaussian  shape 
horizontal  beamwidth.  This  is  a  correction  based  on  Probert- Jones 
correction  for  a  conical  or  elliptical  beam  in  the  azimuthal  beam 
only . 


The  values  in  Table  3  give  a 
stated  above)  of  3.52  x  10^.  The  J^G^  A4)  = 
fore,  53g^  H  =  2.0  2  x  10^^. 


(with  the 
5.75  X  10^. 


loss 

There- 


Zgf  f  (mm^/m^) 


=  1.7  X  10 


(4.2) 


Where  range  is  in  nmi  and  power  received  is  in  milliwatts. 

During  the  Oklahoma  experiment  the  ARSR- ID  was  posi¬ 
tioned  with  maximum  gain  at  3/4®.  This  positioning  directs  some 
energy  onto  the  ground.  Therefore,  ground  reflections' may  be 
present  in  the  ARSR- ID  power  received.  These  ground  reflections 
were  studied  and  will  be  discussed  later.  However,  in  the  compu¬ 
tation  of  presented  herein,  no  attempt  was  made  to  correct 

for  ground  reflections.  In  general  practice,  for  meteorological 
radars,  calibration  and  reflectivity  threshold  levels  are  set 
without  regard  to  antenna  tilt.  Therefore,  it  was  decided  to 
keep  the  computation  of  ^ ^ff  close  as  possible  to  procedures 
used  in  common  practice  by  operating  meteorological  radars. 

The  radar  equati<jn  used  for  the  WSR-57  is  the  equation 
generally  used  in  radar  meteorology.  As  was  discussed  in  the 
section  on  simulation,  the  WSR-57  can  underestimate  the  maximum 
in  a  storm  due  to  the  antenna  beamwidth  and  tilt;  however,  for 
sake  of  consistency  with  the  radar  meteorology  community,  we 
have  analyzed  the  data  using  the  conventional  equation.  The 
radar  equation  for  the  WSR  57  is: 


Pari.meter 


ARSR- ID 


WSR-57 


Peak  transmitted  power  -  P^ 

4700  kW* 

370  kW 

Range  resolution  -  ct/2 

300  meters 

600  meters 

Maximum  gain 

34.3  dB 

38.1  dB* 

Integrated  gain  - 

(integrated  from  -3  dB  point) 

57.6  dB 

60.6  dB 

Horizontal  beatiwidth 

1.35° 

2.0°* 

Vertical  beamwi.dth 
(half  power) 

~6° 

(shaped  to  44°) 

o 

o 

Frequency 

1335  MHz 

2875  MHz* 

Function  of  index  of  refraction 
of  water  - 

.93 

.93 

Correction  for  log  averaging 

2.54  dB 

2.54  dB 

Probert- Jones  correction  for 
Gaussian  shaped  antenna 

0.7  dB 

(azimuth  only) 

1.4  dB 
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24  2 

f.  ,  3.43  X  10  R  P 

ZgCnrniVm'^)  =  - 1^- - - 


(4.3) 


^SR 


t  I’'!"  »♦  S 


(2  )ln  2)  256  X‘ 


where  P^  and  p^  are  in  watts;  6  and  (()  are  in  radians;  t  is  in 
seconds;  X  is  in  meters;  and  R  is  in  nautical  miles. 


and 


Using  the  Parameters  given  in  Table  3,  K^gj^  =  4.7  x  10^' 


Zg  =  7.3  X  lo"^  R^  Pj.  (4.4) 

where  R  is  in  nmi  and  power  received  is  in  milliwatts.  Eqs.  (4.2) 
through  (4.4)  were  used  to  compute  the  results  that  are  discussed 
in  Sections  5  and  6.  Neither  or  K^gj^  contain  atmospheric 

losses  due  to  oxygen  and  water  vapor  because  these  losses  are 
dependent  on  range.  However,  it  should  be  mentioned  that  within 
the  operating  ranges  for  the  storm  cells  analyzed,  the  relative 
difference  between  these  losses  for  the  ARSR  (i.e.,  L-band)  and 
the  WSR  (i.e.,  S-band)  was,  on  the  average,  .3  dB. 

4 . 5  Overlapping  Pulse  Volumes 

The  computer  progreun  and  its  logic  flow  are  discussed 
in  Appendix  A.  However,  since  the  determination  of  those  areas 
for  which  the  two  radars  overlap  in  space  is  important,  we  will 
discuss  it  here. 

The  relative  position  of  the  two  radars  was  shown  in 
Figure  3-1.  Suppose  at  a  point  in  the  program  we  have  averaged 
precipitation  reflectivities  for  both  radars  along  with  ranges 
and  azimuths  to  the  center  of  the  averaged  areas  for  a  selected 
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storm  cell.  We  determine  whether  a  given  ARSR  area  was  overlapped 
by  a  corresponding  WSR-57  area.  To  do  this  we  translate  all 
WSR-57  rangci,  ^^gR'  azimuths,  to  the  ARSR-lD  locations 

to  obtain  Rj'jgp  and  These  R^gp  and  values  are  tested 

to  see  if  they  fall  within  some  small  area  around  the  ARSR-lD 
coordinate  If  and  are  within  this  toler¬ 

ance  area,  the  areas  are  assumed  to  be  overlapping  in  space  (see 
Figure  4-2)  .  The  variables  and  were  inputs  to  the  program 

and  were  chosen  to  give  the  maximum  overlap.  The  amount  of  over¬ 
lap  varies  from  storm  cell  to  storm  cell  and  area  to  area  depending 
on  angular  Jocation  of  the  data  with  respect  to  the  two  radars. 

In  the  cases  where  WSR-57  elevation  scans  were  averaged, 
the  overlap  was  established  for  each  scan  and  the  reflectivities 
for  a  given  area  were  averaged  in  elevation. 


Contouring  and  Area  Calculation 


The  main  IBM  360/91  program  used  for  data  reduction 
printed  the  mean  reflectivities  at  each  scan,  the  rms  reflectivi¬ 
ties,  and  the  average  reflectivities  in  dBZ  for  each  scan  of  the 
WSR  and  ARSR.  The  printout  was  in  B-scan  display  such  as  those 
shown  in  tht;  next  section,  where  range  and  azimuth  are  printed 
in  a  rectangular  grid.  These  B-scans  were  then  stored  on  a 
Hewlett-Packard  9825  calculator  tape,  from  which  a  contoured 
PPI  type  display,  using  NWS  reflectivity  levels  (see  Table  4) , 
was  generattjd.  In  addition  to  the  contours  for  each  set  of 
storm  cell  data,  the  area  with  a  given  reflectivity  level  was 
computed  for  each  storm  cell.  For  example,  consider  a  level  2 
which  has  rc?f lectivi ties  from  30  to  40  dBZ.  All  ARSR  areas 
having  reflectivities  within  this  range  were  summed. 
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Illustration  of  ARSR  and  WSR  overlapping  areas 


^AZ 


ARSR  area 
center 


Figure  4-2.  Illustration  of  overlapping  pulse  volumes 
for  WSR-57  and  ARSR- ID  in  Oklahoma. 
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Z  >  57 
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5.0  DATA  ANALYSIS  -  ENROUTE 


5.1  Contour  Comparison 


Of  the  data  collected  in  Oklahoma,  six  storm  cells  were 
chosen  for  digitization  and  processing.  These  data  were  chosen 
due  to  storm  height,  intensity,  and  geometry.  Other  storm  cells 
were  available  but  they  were  not  digitized  by  the  FAA  because  of 
time  limitations. 


First  the  WSR-57  and  ARSR-lD  data  were  averaged  and 
B-scan  prints  made  to  check  the  quality  of  data  and  geometric 
position  of  the  radar  echoes.  These  B-scans  were  then  calibrated 
and  printed  with  each  radar  in  its  respective  coordinate  system. 
Both  of  these  stages  of  the  data  were  for  checking  the  data  with 
PPI  photographs  for  position  and  reflectivity  levels.  After 
assurance  that  the  data  were  processable,  they  were  run  through 
the  overlap  program  described  previously.  This  program  gives 
B-scans  of  each  radar  in  the  coordinate  system  of  the  ARSR-ID. 

A  B-3can  presents  a  distorted  view  of  a  PPI  section  (range  and 
azimuth)  because  it  presents  polar  coordinates  in  rectangular 
format.  Nonetheless,  it  is  a  simple  method  of  obtaining  radar 
data  in  computer  printout  form.  Figure  5-1  gives  an  example  of 
this  B-scan  presentation.  The  data  in  these  B-scans  were  pro¬ 
cessed  on  an  IBM  360/91.  These  reflectivity  levels  were  then 
put  on  a  cassette  tape  where  they  were  further  processed  to 
give  storm  cell  contours,  areas,  etc. 


Storm  Cell  Comparison  (contours) 


The  B-scan  plots  in  dBZ  were  analyzed  for  the  six  storm 
cells;  processed.  The  data  indicated  that  on  the  average  the 
ARSR--1D  reflectivities  were  approximately  4  dBZ  higher  than  the 
WSR  reflectivities.  No  analytical  justification  for  this  over¬ 
all  difference  in  level  could  be  found,  however,  explanations 
have  been  pursued.  One  possible  source  of  this  difference  could 
be  Cc libration  error  either  in  the  WSR-57  or  the  ARSR-lD.  This 
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error  source  can  possibly  be  ascertained  from  sun  strobe  data 
which  are  available  in  video  tape  form  but  have  not  been  digitized. 
Another  suggested  source  comes  from  the  fact  that  the  maximum  gain 
of  the  antenna  was  positioned  at  3/4®  elevation  during  the  experi¬ 
ment.  It  is  possible  that  this  positioning  could  cause  reflec¬ 
tions  from  the  surface  which  would  increase  the  received  power  of 
the  ARSR-ID.  In  order  to  obtain  some  indication  of  the  magnitude 
of  this  effect,  a  theoretical  analysis  was  performed  with  the 
PMG  of  the  ARSR-lD  at  3/4®.  The  analysis  assumed  that  all  energy 
directed  to  the  ground  would  be  reflected  (i.e.,  a  reflection 
coefficient  =  1) .  Computations  were  performed  at  various  ranges 
and  the  increase  in  total  received  power  due  to  complete  ground 
reflections  determined.  The  results  show  (see  Appendix  D)  that 
for  the  ARSR-ID  the  magnitude  of  the  increased  power  is  depen¬ 
dent  on  range  and  storm  height.  For  storms  with  an  altitude 
extent  greater  than  5  km  (16.4  kft) ,  the  increase  in  power  would 
vary  as  a  function  of  range  from  near  zero  dBm  at  40  km  (21.6  nmi) 
to  near  2.4  dBm  at  200  km  (110  nmi).  A  similar  analysis  for  the 
WSR-57  positioned  at  0®  PMG  would  result  in  an  increase  as  a 
function  of  range  of  approximately  3  dBZ  from  40  km  (21.6)  to 
160  km  (86  nmi) .  A  positioning  of  the  WSR-57  at  or  above  1/2® 

PMG  virtually  eliminates  any  contribution  from  ground  reflections. 

In  practice,  there  is  no  way  to  know  exactly  what  the 
reflection  coefficient  is  for  various  ranges  and  azimuths  at  an 
operating  site.  With  this  in  mind  the  decision  was  made  to  make 
no  adjustments  to  the  experimental  data  due  to  ground  reflections 
for  either  the  ARSR-lD  or  the  WSR-57  although  the  ARSR-ID  was 
positioned  with  the  PMG  at  3/4®  and  the  WSR  was  positioned  at  0® 
for  two  storm  cells. 

As  mentioned  in  Section  4.4,  the  conventional  equation 
for  computing  WSR-57  reflectivity  was  used  for  the  experimental 
data  (i.e.,  maximum  gain  and  2  Zn  2  correction  for  beam  shape). 

Pad  an  integrated  gain  of  the  pattern  for  0®  PMG  been  used  along 
with  complete  ground  reflection  (i.e.,  p  =  1) ,  the  WSR-57 
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reflectivity  values  would  not  have  changed  because  the  3  dB 
increase  from  ground  reflections  is  cancelled  by  a  3  dB  decrease 
from  maximum  gain  to  integrated  gain.  If  the  same  thing  had  been 
done  for  the  ARSR-lD  integrating  the  gain  pattern  from  the  hori¬ 
zon  with  3/4®  PMG,  an  average  decrease  of  1  dBZ  would  result. 

It  appears  that  ground  reflections  cannot  fully  explain 
the  discrepancy  between  the  ARSR-ID  data  and  the  WSR-57  data. 
Despite  this,  the  ARSR-ID  reflectivity  data  presented  has  been 
lowered  by  4  dBZ  to  permit  meaningful  side  by  side  comparisons. 

It  is  felt  that  whatever  the  source  of  this  bias  error,  it 
amounts  to  a  calibration  problem  either  in  the  ARSR-ID  or  the 
WSR-57  and  can  be  accounted  for  in  setting  threshold  levels. 

Thus  it  in  no  way  degrades  the  capability  of  the  ARSR-lD  to 
measure  weather. 

In  Figure  5-2  we  observe  a  storm  cell  comparison 
centered  at  38  nmi  with  respect  to  the  ARSR-ID.  This  storm  cell 
was  recorded  on  September  28,  1977.  The  storm  characteristics, 
as  seen  from  three  elevation  scans  through  the  WSR-57,  indicated 
it  had  maximum  reflectivities  near  the  surface  and  decreased  in 
intensity  with  altitude.  The  contours  in  Figure  5-2  represent 
an  average  of  5  digitized  azimuth  values  and  6  digitized  range 
values.  These  averaged  quantities  due  to  a  digitization  rate  of 
.869  usee  and  the  ARSR-ID  scan  rate  result  in  .5®  azimuth  and 
782  meters  range  resolution.  This  .5®  resolution  in  azimuth  is 
less  than  the  ARSR-lD's  1.35®  and  the  WSR-57  2,0®  azimuthal  beam- 
widths.  This  means  that  the  values  as  presented  do  not  represent 
spatially  independent  estimates.  Averaging  in  azimuth  up  to  the 
horizontal  beamwidth  does  not  substantially  change  the  contours 
but  does  smooth  the  contours  and  reduces  the  area  of  maximum  Z. 
This  is  shown  in  Figure  5-3  for  storm  cell  1  where  both  radars 
have  been  averaged  over  2®  of  azimuth  and  782  meters  range.  The 
general  shape  of  the  cells  is  the  same  in  both  Figures  5-2  and 
5-3.  Storm  cell  1  (Figure  5-2)  shows  a  level  of  6  for  the  ARSR-lD; 
howevei ,  the  WSR-57  does  not  show  this.  The  other  levels  for  this 


Level  Keflectfvity  (dBZ)  Level  Reflectivity  (dBZ) 
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Figure  5-2.  Storm  cell  1  -  ARSR-lD  and  WSR-57  contoured  reflectivity. 
.5®  azimuth  and  782  meter  range  resolution. 
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storm  cell  agree  quite  well.  It  is  possible  that  this  core  was 
at  altitude  and  existed  for  a  very  short  period  of  time,  which 
would  mean  that  the  WSR-57  possibly  missed  it  at  1°  elevation 
angle  and  it  did  not  exist  before  and  after  when  the  WSR  was 
scanning  at  higher  elevations. 

As  mentioned  previously,  all  cells  were  scanned  in 
elevation  (see  Figure  5-4).  Storm  cell  1  had  echoes  up  to  7°. 

The  storm  had  its  highest  reflectivities  at  1®  elevation  with 
the  level  of  the  core  reflectivity  decreasing  about  3  dB  at  3® 
and  decreasing  steadily  from  that  point.  The  last  echoes  seen 
at  7“  elevation  and  39  nmi  indicate  a  storm  height  between  8  and 
10  )cm. 

Storm  cell  2  (see  Figure  5-5)  shows  the  detection  of  a 
cell  at  relatively  long  ranges  (i.e.,  185  )cm  or  100  nmi).  This 
cell  is  interesting  because  echoes  showed  on  the  WSR  at  only  1® 
elevation.  This  means  that  the  storm  did  not  extend  above  9  km 
(29  ):ft)  .  If  we  assume  that  the  storm  extended  to  9  )an,  then 
the  ARSR-ID  with  a  6®  (3  dB)  beamwidth  has  less  than  three  quarters 
of  the  beam  filled.  The  simulation  in  this  case  would  predict 
the  ARSR  to  be  at  least  1.2  dBZ  below  the  WSR.  This  1.2  dBZ  would 
be  difficult  to  see  in  the  contours;  however,  there  is  no  indi¬ 
cation  that  the  ARSR  is  lower.  Instead  the  contours  from  the  WSR 
and  ARSR  are  approximately  equal. 

Storm  cells  3  and  3A  (Figures  5-6  and  5-7)  are  the  same 
storm  separated  in  time  by  about  4  minutes.  Again  the  general 
shapes  of  echoes  are  similar.  The  ARSR  area  at  level  2  seems  to 
be  smaller.  It  is  judged  that  the  information  as  presented  by 
the  ARSR- ID  would  be  as  helpful  to  air  traffic  controllers  as 
that  given  by  the  WSR. 

Figure  5-8  shows  a  rather  large  and  complex  storm  cell. 

Again  the  similarities  are  apparent.  This  cell  can  be  used  to 
illuitrate  the  possiblity  for  implementation  of  a  2  level  display 
by  the  FAA  for  controller  use.  Figure  5-9  illustrates  where  the 


Storm  cell  3A  -  ARSR— ID  and  WSR-57  contoured  reflectivity 
.5®  azimuth  and  782  meter  range  resolution. 


Reflectivity  (Z) 
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Example  of  contoured  reflectivity  showing  two  levels. 
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thresholds  could  be  set  for  an  installation;  30  dBZ  (level  2)  and 
41  dBZ  (level  3  and  above) .  The  information  shown  by  the  two 
radars  would  be  virtually  the  same. 

A  stratiform  rain  condition  is  not,  in  general,  hazardous 
to  aircraft.  However,  storm  cell  5  is  shown  (Figure  5-10)  to 
indicate  how  this  type  rain  might  appear  on  the  scope.  There  are 
no  reflectivities  above  level  2  and  the  display  appears  patchy. 

In  this  particular  rain,  which  actually  is  not  an  intense  storm, 
there  were  no  reflectivities  above  37  dBZ.  In  general,  this  low 
level  patchy  signature  is  indicative  of  stratiform  precipitation 
and  is  not  difficult  to  recognize  after  seeing  it  several  times. 

5 . 3  Area  Comparisons 

To  yield  a  measure  of  the  comparability  of  the  results 
from  two  radars,  area  comparisons  were  calculated.  For  a  given 
storm  cell,  the  area  included  at  each  level  was  computed  for  the 
two  radars  for  all  six  cells.  The  ratios  of  the  ARSR-ID  area  to 
the  WSR  area  are  shown  in  Figure  5-11.  The  level  1  indicated  in 
this  figure  is  different  from  the  NWS  level  1  in  that  it  indicates 
dBZ  values  between  20  and  30  dBZ,  There  is  quite  a  bit  of  scatter 
in  the  data.  However,  the  mean  is  very  close  to  1.0.  This  result 
supports  the  comparability  of  the  results. 

5 . 4  Comparison  with  the  Simulation 

In  theory  the  simulation  indicates  that,  given  the 
storm  profile  at  some  range,  one  can  predict  what  the  ARSR-ID 
should  measure.  We  have  chosen  18  profiles  from  the  WSR  measure¬ 
ments  in  several  storms.  This  was  done  by  talcing  an  average 
reflectivity  within  a  .7®  by  1  km  area  at  each  WSR  elevation 
angle.  This  reflectivity  was  then  assigned  to  an  altitude  cor¬ 
responding  to  the  mid-point  of  the  WSR  beam  at  that  range.  It 
should  be  noted  that  these  profiles  obtained  from  the  WSR- 57  are 
averages  of  the  true  profile.  That  is,  the  scatterers  within  the 
WSR-57  2.0®  beam  have  been  averaged.  Technically,  the  ARSR-ID 
averaged  the  true  profile  over  its  beamwidth.  Therefore,  since 


Level  number 


Figure 


Storm  intensity  levels 


-11.  Ratio  of  ARSR-lD  area  to  WSR-57  area  at  a 
level  vs  level  number.  The  solid  line  is 
the  average  of  all  points. 
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the  "true  profile"  was  not  put  into  the  simulation,  we  should  not 
expect  the  ARSR-lD  measured  reflectivity  to  be  exactly  equal  to 
the  simulation.  However,  the  two  reflectivities  should  be  close  if 
the  simulation  is  representative  of  the  actual  ARSR-ID  measurement. 
Figure  5-12  presents  the  results  in  a  regression  plot.  Each  pro¬ 
file  was  run  through  the  simulation  program  and  an  ARSR  reflectivity 
computed.  There  is  a  .89  correlation  indicated.  Th linear 
regression  equation  gives  a  slope  of  .9  between  the  measured  and 
the  simulated  ARSR  values.  Perfect  agreement  would  be  a  slope 
of  1  and  intercept  of  0 . 

The  ARSR- ID  measured  values  were  also  compared  to  the 
maximum  V7SR  measured  reflectivities  for  the  18  profiles.  These 
maxima  were,  in  general,  at  the  lowest  elevation  angle  0®  or  1® 
for  the  V/SR-57.  Figure  5-13  gives  these  results.  The  slope  of 
.97  and  intercept  of  -.88  from  the  regression  are  close  to  inter¬ 
cept  0  and  slope  1  which  perfect  agreement  would  give. 

5.5  Comparison  of  MTI  and  Log  Video  Data 

All  reflectivity  values  previously  discussed  have  been 
computed  from  log  video.  During  the  Oklahoma  experiment  Moving 
Target  Indicator  (MTI)  data  were  also  recorded.  The  MTI  mode  is 
used  on  FAA  radars  to  eliminate  fixed  targets  such  as  ground  clut¬ 
ter.  To  accomplish  this  the  velocity  notch  of  the  MTI  filter  is 
centered  at  zero  velocity.  As  a  result  slowly  moving  weather 
targets  may  also  be  eliminated.  However,  in  many  thunderstorms 
internal  velocities  may  be  sufficiently  high  to  fall  in  the  pass 
band  of  the  filter  and  allow  the  weather  returns  to  be  detected. 

For  this  reason  MTI  video  data  were  recorded  simultaneously  with 
log  video  data.  The  resulting  data  showed  that  portions  of 
storms  were  detected  on  MTI,  however,  the  reflectivity  measure¬ 
ments  were  not  good  for  two  reasons.  On  the  edge  of  the  storm 
cells,  where  reflectivities  were  below  30  dBZ,  there  were  differ¬ 
ences  between  MTI  and  log  ohat  were  inconsistent.  That  is,  the 
MTI  at  tiiiies  measured  patt<;rns  and  reflectivities  similar  to  the 
log  video  and  other  times  they  were  dissimilar.  In  those  areas 
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ARSR-1D  theoretical  reflectivity  (dBZ) 


Figure  5-12.  Measured  ARSP-lD  reflectivity  vs  simulated 
ARSR-lD  reflectivity  for  18  WSR  measured 
profiles . 


ARSR-1D  measured  reflectivity  (dBZ) 
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of  reflectivity  greater  than  approximately  30  dBZ,  limiting  in 
the  MTI  receiver  caused  the  signals  to  saturate  so  that  the 
maximum  reflectivities  could  not  be  accurately  measured.  Figure 
5-14  shows  an  example  of  how  the  typical  MTI  to  log  video  com¬ 
parison  appeared  where  the  solid  symbols  indicate  areas  of 
saturation.  This  particular  cell  was  one  that  showed  a  good 
pattern  at  level  2  reflectivity,  but  reflectivities  greater  than 
level  3  could  not  be  accurately  measured. 

The  results  indicate  that  the  MTI  of  the  ARSR  radar 
as  presently  designed  is  not  good  for  weather  reflectivity 
measurements,  due  to  limiting  of  the  receiver  and  the  zero  vel¬ 
ocity  notch  of  the  filter.  The  MTI  weather  reflectivity  data 
analyzed  was  beyond  ground  clutter  ranges.  Over  ground  clutter 
MTI  weather  reflectivities  would  be  further  suppressed.  Further 
work  may  be  desirable  for  determining  how  to  measure  MTI  weather 
reflectivities . 


6.0  DATA  REDUCTION  AND  ANALYSIS  -  TERMINAL 

Considerable  effort  was  put  into  the  data  reduction 
phase  of  data  taken  in  the  New  Orleans  experiment.  However, 
proolems  were  encountered  which  have  not  been  corrected  to  this 
date.  These  problems  were  from  numerous  sources.  In  some  cases 
calibrations  were  in  error  and  in  others  there  were  spikes  in 
the  raw  data  that  could  not  be  removed.  All  problems  pointed  to 
difficulties  in  the  recording  of  the  data  and  experimental  pro¬ 
cedures.  In  no  case  was  there  an  indication  that  the  problems 
stemmed  from  the  radars.  Due  to  these  problems,  no  data  from  the 
terminal  experiment  are  available  for  inclusion  in  this  report. 
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7 . 0  INTERPRETATION  OF  RESULTS 
7 . 1  Enroute  Radars 

The  data  analyzed  from  the  Oklahoma  experiment  and  the 
simulation  have  indicated  that  the  FAA  enroute  radars  are  capable 
of  measuring  reflectivity  sufficiently  accurately  to  provide  a  con 
toured  weather  display  for  air  traffic  controllers.  The  experi¬ 
mental  results  presented  herein  have  a  bias  of  app  'o;  .^mately  4  dBZ 
which  has  not  been  accounted  for  to  date.  Since  this  is  a  bias 
error,  the  ARSR  radars,  provided  they  are  properly  calibrated, 
should  measure  storm  intensity  and  position  as  well  as  the  NWS 
WSR-57,  at  ranges  within  185  km  (100  nmi) . 

The  simulation  indicates  that  with  a  3/4®  PMG  the 
ARSR- ID  should  detect  as  well  as  the  WSR  for  convective  storms 
when  detecting  within  240  km  (130  nmi),.  because  within  these 
ranges,  heights  for  severe  storms  should  be  sufficient  to  fill 
most  of  the  beam.  Beyond  this  range  some  statistical  correction 
for  beamfilling  and  earth  curvature  effects  could  be  applied  to 
the  received  signal  to  estimate  the  true  reflectivity.  The 
amount  of  correction  will  depend  on  the  radar  and  will  be  more 
for  the  ARSR- ID  with  its  6°  (3  dB)  beamwidth,  than  for  the  ARSR-2 
and  ARSR-3  which  have  a  4®  beamwidth.  The  WSR-57  will  also 
measure  signals  that  are  degraded  by  beamfilling  and  earth  curva¬ 
ture  effects  if  operated  beyond  240  km.  A  statistical  correction 
should  also  be  applied. 

It  should  be  emphasized  that  the  experiment  in  Oklahoma 

was  conducted  with  the  ARSR- ID  operated  in  modes  suitable  for 

weather  detection.  That  is,  linear  polariza '-.i on  was  used  in  all 

six  cell  storms  analyzed  here.  Where  the  STC  curve  used  was 
-  4 

equivalent  to  R  ,  an  adjustment  was  made  in  the  computer  program 

_  2 

to  correct  back  to  an  R  STC  curve.  This  indicates  that  th-: 
enroute  radars  must  be  operated  in  an  optimum  weather  mode  in 
order  to  measure  reflectivity  accurately.  This  may  require  modi¬ 
fications  and/or  adjustment  to  the  radar  or  correction  to  the 
received  weather  data  to  produce  representative  reflectivity 
estimates . 
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7 , 2  Terminal  Radar 

The  simulation  has  indicated  that  the  terminal  radar 
ASR-8  low  beam  and  ASR-5  should  have  no  problem  detecting  accurate 
reflectivity  data  when  operated  within  conventional  60  nmi  range. 
Although  clutter  was  not  addressed  in  the  simulation,  the  basic 
assumption  is  that  weather  measurements  inside  of  clutter  ranges 
are  not  accurate. 
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8 . 0  CONCLUSIONS 

The  following  conclusions  are  made  based  on  the  results 
from  this  study. 

(1)  All  ARSR  radars  can  be  used  to  measure  and  display 
meaningful  weather  data  when  operated  in  a  mode  optimized  for 
weather  detection  within  240  km.  This  conclusion  is  based  on 
both  experimental  results  from  the  ARSR- ID  and  simu^tuod  results. 
The  240  ]:m  range  restriction  is  imposed  because  beyond  this 
range,  the  effects  of  earth  curvature,  partially  filled  beams, 
and  intecjration  greatly  degrade  the  estimates  of  storm  intensity. 

(2)  All  ARSR  radars  and  the  WSR-57  will,  in  general, 
underestimate  the  maximum  reflectivity  with  the  amount  of  under- 
estimatii^n  dependent  on  storm  range  and  position  of  the  maximum 
in  altitude.  For  storms  with  maximum  reflectivity  on  the  ground 
the  underestimation  increases  with  increasing  3  dB  beamwidth. 
Thus,  the  WSR-57  with  the  2  degree  beam  should  be  closer  to  the 
maximum.  When  the  maximum  is  positioned  at  high  altitudes  the 
ARSR  radars  should  detect  the  storm  intensity  more  accurately 
than  the  WSR-57  radar  in  a  normal  scan  mode  (i.e.,  .5®  or  1® 
elevation  angle) .  This  result  is  significant  since  there  is 
evidence  that  severe  storms  have  maximum  intensity  at  higher 
altitudes  during  the  period  of  peak  development. 

(3)  The  ARSR  and  the  WSR-57  radars  measure  a  reflectivity 
closer  to  an  average  of  the  total  profile  reflectivity  (i.e., 

<Z>)  than  to  the  maximum. 

(4)  The  ARSR  MTI  systems  as  they  are  presently  configured 
arc  not  .suitable  for  accurate  weather  measurement. 

(5)  At  ranges  greater  than  240  Ion  the  ARSR  radars  should 
have  a  statistical  correction  factor  applied  to  the  received 
signa]  tC'  account  loi  partial  beamfilling,  earth  curvature  and 
integration  effects.  It  siiould  be  noted  that  the  WSR-57  should 
also  have  such  correction  if  operated  beyond  this  range.  The 
exact  nature  of  this  statistical  correction  would  have  to  be 
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determined  based  on  a  study  of  average  heights  and  profiles  of 
severe  storms.  Since  such  a  correction  factor  would  be  statis¬ 
tical  in  nature,  it  would  not  eliminate  all  reflectivity 
measurement  errors  at  long  ranges  but  rather  reduce  them. 

(6)  The  simulation  indicates  that  the  ASR-8  low  beam  and 

ASR-5  are  capable  of  detecting  weather  accurately  within  their 

normal  operating  range  and  outside  the  limits  of  ground  clutter 

(i.e.,  MTI  area) .  This  statement  assumes  an  operating  mode  which 

is  suitable  for  detecting  weather  (i.e.,  linear  polarization  and 
-2 

R  STC  curve)  . 

(7)  The  ASR-8  high  beam  is  not  suitable  for  weather 
detection  when  operated  as  the  sole  source  of  the  measurement. 
However,  when  it  is  used  in  a  range  gating  mode  with  the  low 
beam  and  restricted  to  close-in  ranges,  it  may  measure  reflec¬ 
tivity  as  accurately  as  the  low  beam.  The  reason  for  restriction 
of  the  high  beam  to  short  ranges  is  that  its  maximum  gain  is 
positioned  at  a  higher  elevation  angle  which  causes  it  to 
illuminate  the  upper  portions  of  storms  at  shorter  ranges  than 
the  low  beam.  The  high  beam  is  operationally  utilized  at  close- 
in  ranges  to  minimize  the  reception  of  ground  clutter.  This 
improves  MTI  and  weather  reflectivity  over  ground  clutter 
performance , 

(8)  The  location  and  strength  of  reflectivity  echoes 
should  aid  air  traffic  controllers  in  determining  areas  of 
intense  precipitation.  It  should  be  noted,  however,  that  research 
studies  concerning  turbulence  associated  with  thunderstorms  has 
not  established  a  clear  relationship  between  intense  reflectivity 
areas  and  turbulent  areas.  In  fact,  areas  of  significant  turbu¬ 
lence  have  been  observed  to  exist  in  the  clear  air  outside  the 
areas  of  reflectivity  associated  with  the  thunderstorm.  Based 

on  such  studies,  the  recommended  criteria  is  for  aircraft  to 
avoid  higher  reflectivity  (46  dBZ  and  above-levels  4,  5,  and  6) 
storm  cores  by  10  to  15  miles.  Also,  it  can  be  said  that  high 
reflectivity  areas  (50  dBZ  and  above-levels  5  and  6)  will  contain 
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intense  precipitation  and  often  aircraft  damaging  hail. 
Therefore,  if  such  reflectivity  data  is  properly  interpretated, 
such  displayed  information  could  materially  assist  the  air 
traffic  controllers. 

NOTE:  Varying  amounts  of  equipment  modifications  and/or 

additions  would  be  required  (depending  on  the  Model  FAA  radar) 
for  simultaneously  obtaining  optimized  displays  of  aircraft 
and  calibrated  weather. 
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THE  COMPUTER  PROGRAM 


1.0  INTRODUCTION 

For  the  study  of  severe  weather  detection,  an 
experiment  was  conducted  in  Norman,  Oklahoma,  during  part  of 
September  and  October  of  1977,  as  described  in  the  body  of  this 
report.  Weather  data  were  collected  by  the  ARSR  radar  and  the 
WSR-57  radar  simultaneously.  In  order  to  compare  the  ARSR  data 
with  the  WSR-57  radar,  an  IBM  360/91  computer  program  was  written 
to  reduce  the  data  from  both  radars.  A  description  of  this 
computer  program  is  given  in  this  Appendix  and  the  associated 
flow  chart  for  the  program  is  given  in  Table  A-9 .  The  program 
described  herein  is  for  the  enroute  data  reduction.  The  program 
for  terminal  data  reduction  is  similar  with  only  slight 
modifications . 

2.0  DESCRIPTION  OF  THE  PROGRAM 

In  order  to  compare  the  radar  reflectivity  factors  of 
the  ARSR  and  those  of  the  WSR-57,  overlapping  areas  of  radar 
video  are  digitized  and  recorded  on  VQR  tape  by  NAFEC.  The  com¬ 
puter  program  described  here  performs  the  following  functions. 
First,  it  reads  the  ARSR  data  from  the  VQR  digital  tape  and 
obtains  a  data  matrix  of  reflectivity  in  dBZ.  It  follows  by 
reading  the  WSR-57  data  from  another  VQR  tape  and  generates  a 
data  matrix  of  reflectivity  in  dBZ.  The  overlapping  area 
between  the  two  radars  is  selected  and  a  data  matrix  of  WSR-57 
reflectivity  values  with  respect  to  the  ARSR  coordinate  system 
is  obtained. 

2 . 1  Read  ARSR  Data 

The  main  progreuti  first  reads  in  all  the  necessary  input 
parameters  for  the  ARSR  data  on  IBM  cards.  A  brief  description 
of  these  parameters  is  as  follows:  (1)  The  ARSR  calibration 
curve  of  tape  units  (i.e.,  digitized  voltages)  versus  power 
received  (dBm)  with  noise  removed.  The  reason  for  removing  the 
noise  is  as  follows.  Due  to  the  digitization  process,  the  noise 
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of  the  data  in  terms  of  tape  units  on  the  VQR  tape  is  different 

from  "he  noise  in  tape  units  of  the  calibration.  In  order  to 

use  t'le  calibration  curve,  the  data  and  the  calibration  must 

have  the  same  reference  for  noise.  This  is  accomplished  by 

removing  the  respective  noise  in  tape  units  (TU)  from  both,  so 

the  nc ise  is  at  0  TU  for  both  the  data  and  the  calibration. 

The  actual  procedure  for  the  computation  is  explained  in  a  later 

section.  (2)  The  STC  curve  attenuation  in  dB  versus  range  (nmi) 

-2  -4 

is  the  next  input  item.  Either  the  R  or  R  curve  is  used. 

_  2 

In  Figure  A-1  note  the  R  curve,  the  attenuation  starts  from 

32  dB  at  range  1  nmi  and  decreases  to  1  dB  at  80  nmi.  And  for 
-4 

the  R  curve,  the  attenuation  starts  from  60  dB  at  1  nmi  and 
decreases  to  1  dB  at  50  nmi.  The  input  also  includes  K,  the 
constant  used  in  the  radar  equation,  and  VNOISE,  the  noise  level 
(tape  units)  of  the  data.  The  neighborhood  in  terms  of  azimuth 
in  degrees,  range  in  nmi,  inputs  used  to 

check  for  overlapping  areas.  Finally,  NO,  the  number  of  range 
bins,  and  NAR,  the  number  of  azimuth  sweeps  to  be  averaged,  are 
also  read  in.  In  those  cases  where  range  exceeds  approximately 
60  nmi,  a  range  correction,  RCORR,  is  read  in  to  obtain  the  true 
range  of  the  window.  This  step  is  required  because  the  VQR 
machine  can  only  digitize  60  nmi  of  range  extent.  There  is  a 
difference  between  the  true  north  and  the  AGP  (azimuth  change 
pulse)  which  was  used  in  recording  the  video  during  the  experi¬ 
ment.  Therefore,  this  azimuth  correction,  AZ-REL,  is  read  in  to 
generate  the  true  azimuth  for  the  VQR  window.  The  format  and  an 
example  of  these  inputs  are  given  in  Table  A-1.  After  all  the 
above  parameters  are  read  in,  the  main  program  calls  the  sub¬ 
routine  WSSCAN  to  read  in  the  digital  data. 


Before  proceeding  to  the  subroutine  WSSCAN,  we  give  a 
brief  description  of  how  the  data  are  stored  on  tape.  The  data 
on  the  VQR  tape  are  from  several  windows  or  scans  which  are 
separafed  by  end  of  file  marks.  The  VQR  tape  format  is  given  in 
Table  \-2.  Each  scan  is  made  up  of  a  number  of  VQR  azimuth  sweeps 
where  each  sweep  covers  from  R^  to  R2  in  range.  In  front  of  each 


Card  ^  j. 

Input 

ARSR 

WSR-57 

1 

Calibration  of  tape  units 
vs  power  received  (dBm) 

1 

0 

L, 

File  No. 

1.1 

2.2 

3 

NO 

6 

6 

4 

NAR 

5 

5 

5 

AAZ 

(azimuth  sweep  resolution) 

0.1 

0.12 

6 

AZ-REL 

(azimuth  correction) 

13.74° 

0.0 

7 

RCORR 

(range  correction) 

0.0  nmi 

57.44 

8 

^AZ 

0.3° 

.3° 

9 

0 . 5  nmi 

.  5  nmi 

10 

S 

1.67E16 

4.5E16 

11 

VNOISE 

90  (TU) 

233 

12(a~) 

STC  (R"^)  curve 

dB  vs  range  nmi 

None 

13(a~^ 

STC  (R  ^)  curve 

Same 

None 

14 

channel 

CHANL  =  ’ 2 ’  ; 
or  '1' 

CHANL  =  ' 2 ' ; 
or  '1' 
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Preamble ; 


Table  A-2  (.continued) 
VQR  TAPll  FORMAT 

word  1  Tape  identity  ^digital) 


6  lower  half  word  contains  the  scan  number 

7  command  word  0 

8  command  word  1 


command  word  0 


29 

28  27 

26  25 

24 

23 

15 

14 

1 

0 

0 

SI 

c 

no ,  of 

computer 

PD 

words  1 

jer  sweep 

packing  density  is  3-10  bit  character 
SI  sample  interval  0.869  usee 

\  =  no.  of  com])uter  words  per  sweep 


command  word  1 


12  11 


r;uige  start 


azimutli  start 


used  to  chock  tor  preamble 


Page  Seven-A 


rable  /\-2  (continued) 
VQR  TAPP  I'OIMAT 


Viceo : 


computer  word  1 


computer  word  2 


29 

20 

19 

10 

9 

1 

0 

channel 

1 

channel 

2 

channel  1 

sample 

2 

sample 

1 

29 

20 

19 

10 

9 

0 

channel 

2 

channel 

1 

channel  2 

n 

sample  3 

sainjile  2 

bits 


bits 


50  bit 


, start  range  R 


one 

azimuth 
swe  jp 


- AAAAAAV' - - 

sample  M 

sample  2 

- MA/WW' - ^ - 

sample  1 

channel  1 


sample  M 


-AW\/\MA- 


-'v'vAM  AAA- 


sample  2 


sample  1 


channel  2 


1^.869  I 

I''  ysec  n 


where  M  -  tot  il  no.  of  range  bins 

M  -  N  no.  of  computer  words  x  3/2 

...  ...  4024  X  3^ 

AS  -  azimuth  sweep  =  ( - ^ - )/N 

14  =  Rj  +  (M-I3  X  0.809  usec/12.563  usee  per  nmi 


L.icl  s.iiiiple  I  i  video  is  a  digital  number  between  0  and  1024. 
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scan  thei‘"  i  a  ^a'eainVile  which  contains  the  ;;can  nun', her,  start 
range,  start  azin'.uth,  and  the  number  of  computer  v.'ords  er  sweep 
from  v;hich  the  number  of  range  bins  is  ceterrair.ed.  Each  range 
bin  is  0.869  :iSec  long  and  the  azimuth  resolution  is  0.1°  for 
ARSR  data  and  0.12°  for  WSR  data.  The  resolution  cf  the  data  on 
tape  is  too  small  for  comparison  of  the  mean  reflectivities  from 
the  two  radars,  therefore,  the  program  averages  h'O,  number  of 
range  bins,  and  MAR,  number  of  sweeps  (or  azimuths) .  WSSCAN  is 
the  subroutine  which  reads  the  data  from  the  digital  tape  and 
perform:-  the  averaging  properly.  It  returns  a  data  matrix  of 
averaged  tape  units  to  the  main  program. 

The  computation  from  tape  units  to  dBZ  is  performed  in 
the  main  program  in  the  following  manner.  Consider  a  tape  unit 
from  the  data  matrix,  the  VNOISE  in  tape  units  is  subtracted  from 
that  tape  unit,  then  the  power  received  ■' n  dBm  corresponding  to 
that  tape  unit  is  obtained  by  interpolation  on  the  input  calibra¬ 
tion  curve.  For  example,  assume  that  the  noise  level  of  the  data 
and  < he  calibration  are  88  and  115  tape  units,  respectively. 

By  si  btracting  88  TU  from  every  point  in  the  data  and  115  TU  from 
ever'^  point  in  the  calibration  curve,  we  bring  the  noise  level 
of  be  th  parameters  down  to  zero.  Then  the  power  received  for  the 
data  is  obtained  by  applying  the  proper  calibration  value,  which 
is  determined  by  interpolation  from  the  calibration  curve  with 
noise  already  remo-ed.  This  power  received  (dBm)  is  increased 
by  th  amount  of  attenuation  in  d.B  by  interpolating  the  STC  curve 
at  tie  appropriate  range.  Using  this  value  of  power  received, 
the  corresp'ondi ng  dBZ  value  is  obtained  by  using  the  radar 
equation.  The  window  of  ARSR  dBZ  values  is  then  constructed. 

An  example  of  this  data  matrix  is  given  in  Table  ?.-2  .  The  range 
and  t  zimuth  values  corresponds  to  cell  center.  Notice  when  the 
comfM  t  ed  dnz  is  below  20,  the  value  11  is  printed.  When  t.he  data 
is  ai  saturation,  the  v,ilue  99  is  printed  in  its  place. 

2  .  1  F-ead  W-UR-f?  Data 

"h.o  next  sto'  is  t read  in  ail  the  i.nrut  j  arameters 
for  f  lie  Wdjf.-'''’  data.  i’i’ey  are  similar  to  those  foi'  the  AREP  data. 
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The  subroutine  WSSCAN  again  reads  the  WSR  data  in  and  averages 
properly  and  returns  with  a  data  matrix  of  tape  units.  The  VNOISF 
in  tape  units  is  removed  from  each  of  the  elements  of  the  data 
matrix.  By  using  the  proper  WSR  calibraion  curve  and  radar 
equation,  v;e  obtain  the  window  of  WSR-57  dBZ  values.  In  the 
following  section,  the  method  of  obtaining  the  overlapping  area 
of  the  two  radars  is  described. 

2 . 3  Find  Overlapping  Area  of  ARSR  and  WSR  Windows 

Each  of  the  windows  of  data  is  given  with  respect  to  its 
own  radar  in  terms  of  range  and  azimuth.  See  Figure  A-2  for  the 
sketch  of  the  ARSR  and  WSR  windows.  In  order  to  compare  the  mean 
radar  reflectivity  of  the  ARSR  and  WSR-57  using  the  ARSR  radar 

position  as  reference,  the  coordinates  (R  ,  AZ  )  of  the  center  of 

Co  (0 

a  WSR-57  cell  must  be  transformed  into  (R',  AZ')  the  new  coordinates 
with  respect  to  the  ARSR  radar.  The  equations  used  for  the  con¬ 
version  of  coordinates  are  given  in  Figure  A-3.  The  new  corrdinates 
(R',  AZ ' )  are  then  used  to  check  of  the  position  in  question  fell 
within  the  neighborhood  ^Range^  ARSR  cell  indicated 

by  does,  the  radar  reflectivity  of  that  WSR  cell 

is  placed  in  the  position  of  the  data  matrix  of  that  ARSR  cell. 
Hence,  a  data  matrix  with  ARSR  coordinates  is  filled  by  a  scan  of 
overlapping  WSR  reflectivity  values.  An  example  of  this  is  given 
in  Table  A-4.  When  there  is  no  overlap,  zero  is  placed  in  the 
data  matrix.  Table  A-5  gives  the  rms  values  within  an  averaged 
area  for  a  single  scan. 

2 . 4  Example  of  Printout 

The  program  continues  by  reading  the  next  elevation  scan 
of  the  WSR-57  data.  After  printing  out  the  data  matrix  of  that 
WSR-57  scan,  the  average  of  the  mean  WSR  radar  reflectivity  for 
these  scans  and  previous  scans  is  computed  and  the  average  of  the 
rms  of  these  scans  is  also  computed  (see  Tables  A-6  and  A-7) .  As 
each  scan  is  read  in  and  processed,  the  average  of  the  mean  and 
the  rms  of  the  radar  reflectivity  factors  of  the  cumulative  scans 
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are  printed  out.  This  continues  until  all  the  desired  WSR  eleva¬ 
tion  scans  have  been  presented.  Finally,  a  data  matrix  of  the 
differences  between  the  ARSR  dBZ  values  and  the  WSR-57  dBZ  values 
is  printed.  An  example  of  this  is  shown  in  Table  A-8.  This 
completes  the  description  of  the  progreun. 

3.0  THE  FLOW  CHART 

The  flow  chart  of  the  progreun  is  given  in  Table  A-9. 

It  is  divided  into  four  parts:  Read  ARSR  data.  Read  WSR-57  data, 
Find  overlap  area  of  WSR-57,  and  ARSR  data  and  printout  data 
matrices.  The  details  of  the  flow  chart  follow  the  description 
of  the  program  very  closely.  Finally,  the  listing  of  the  program 
is  given  at  the  end  of  this  Appendix. 


EXAMPLE  OF  B-SCAN  SHOWING  DIFFERENCE  BETWEEN  FAA  RADAR  AND  WSR-57  RADAR 
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II  — 


>  n.  <\:  if  O  4 

I  I  ^ 


^  ^  (Sj  ^  c:  i 

I  —  4- 


h-  iT  ^  ^  o  r**- 

•  •  •  •  — • 


c  <V  ^  IT  r- 

•I  •  •  ^ 


tf  r-  if  ^  ^  - 

I 


<n  <r  c  ^  rv.  r*  < 

^11  • 

I 

rv  A  <v  o  ( 

till  —  , 


c  —  r*.  ir> 

—  I  I 


A  ^  ^  'A  c  a: 


«-•  tf  r*  A.  a 

I  I 


^  ^  A  A.  C 

II  It 


G.  <$  \f  r  ^ 

A  A  ^  I 


^  A  IT 

^  ^  I 


^  a  r  N>  c 


a  If  «»  N»  >4 
•-  A  ^  I 


^  4*  4-  — 

t  ^ 

I 


tr  if  ^  ■ 

I  I  —  —  ^  . 

I 

I  -O  A  —  «  ^  • 

I  I 


^  A  '*>  4  r-  -I  c.'  *-«  IT  or  rsi  u^  0 


A  r.#  >  g  M  m  4  t 

fA  1*1  m  i»>  p ,  f\-,  , 


*>  e  .  Ai  r'i  ('I  (A  fA  m  f•^  r^  «a  rA 


3‘;"US>: 
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Table  A- 9 

l-LOW  CilART  OF  THL  COMl'UTER  PROGRAM  TO  PROCESS  1-AA  VQR  DATA 


(1)  Read  window  of 
ARSR  data.  Read 
Input  parameters 


Read  ARSR  calibration  of  TU  vs  power  received 

I 

Read  in  ARSR  file  no. 

Read  NO  -  no.  of  range  to  be  averaged 

NAR  -  no.  of  azimuth  sweep  to  be  averaged 

i 

Read  in  neighborhood  e  (“)  5  e_(nmi) 

I  A  2  K 


Read  in  K  ,  V  (TU) 

p  noise 


Read  in  STC  curve  dB  vs  range  (nmi) 

SUBROUTINE 

WSSCAN  ^ . -  . > 

,  y  \/  from  tape 

Read  ARSR  video  from  VQR  tape 

'I' 

Read  preamble  obtain  scan  no.  §  computer  words/sweep 
Average  tape  units  for  NO  of  range  bins 
Average  NAR  of  azimuth  sweeps 
Obtain  data  matrix  cr  tape  units 


Return  to  main 
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Table  A-9  (continued) 


Generate  azimuth  ii  range  arrav  for  ARSK  data 

i 

Obtain  dBm  from  calibration  curve 

Correct  power  received  by  STC  curve 

Compute  dBZ  from  Radar  Equation 

'1/ 

Print  out  ARSR  dBZ  matrix  of  the  window 

^  (2)  Read  WSR-57 

Read  in  WELE  -  no.  of  elevation  data  window 

I 


I  =  1 


Initialize  arrays  WA,  VAP,  SIG  5  WAP 

'i'  (Read  WSR  input 

Read  in  WSR  file  no.  parameters) 

Read  in  NO  -  no.  of  range  bins  to  be  averaged 

NAR  -  no.  of  azimuth  sweeps  to  be  averaged 

Read  iii  K  and  VIMOISE  for  WSR 

Read  in  WSR  calibration  curve  of  FU  vs  power  received 
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(Read  one  scan  of 
WSR-57  datii  from  Lapo) 


Read  WSR  video  from  VQR  tape 
Obtain  preamble 

'I' 

Average  NO  of  range  bins 
Average  NAR  of  azimuth  sweeps 

'I 

Obtain  WSR  matrix  of  averaged  (TU)  and  matrix  of  RMS  (TU) 


Generate  azimuth  and  range  arrays 

Remove  VNOISE  and  obtain  dBm  from  WSR  calibration 

>1' 

Compute  matrix  of  dBZ  from  averaged  TU 
Compute  matrix  of  RMS  in  dB 

'I 

II  =  ^  (3)  Find  overlap  area 

of  WSR  and  ARSR 
windows 
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SUBROUTINE 
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idPlc  ,\-9  (continued) 


Store  WSR  dBZ  (II,  12)  into  matrix  of  ARSR  grid  of  (13,  14)  position 


SIG  (13,  14)  =  RMS  of  WSR  (II,  12) 

DIFF  (13,  14)  =  ARSR  (13,  14^  -  WSR  (II,  12)  in  dBZ 

12  =  12  -  1 


(4)  Print  out 

data  matrices 


I’rint  uT  ARS  .  i  'i  of  data  matrix  of  WSR  dBZ 


Obtain  and  ■■  int  out  average  of  mean  WSR  dBZ 
vo.  gc  of  the  RMS  of  WSR  scans 


■  •..'l  v"  ■ 
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Table  A  9  (continued) 


1  =  1  +  1 


Obtain  difference  between  ARSR  and  WSR  dBZ  values 
and  print  out  matrix  of  the  difference 


ro'-ip  T  Lt 


COMPUTER  LISTING 

•'T/Ml  :  nPTTnM‘;(vtA  \)! 
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STUPCP  I  T«:TT^'^, 


OP  or 

/oer'i/^r  1  */ 

D'-l  Vi'triTSr  mtm  -  xTFOMftL  : 

nri„  FiArs  rtt  (1)  TMTTfM’P's 

nri  (  TFMfjr^,  TFMoa  ,  ^1  ^AT  HIv: 

nrt  rw(i?n)  C^Aon)  l' TT  m  l  ?om  ^ 

DTL  rwA'^'i  rMap(i)  f  y  t*- pnai,  s 

nri  (  PO  TMTT  (  1?,44)  » 

Tf'TT  (  3?1  .^1  • 

/'7FP  T''JTT  (.71  «  FOFTLnN  MOOM)  TN  A7  •/ 

PfP  TAfT  (0,'i)  /*  Fopri.ON  riOIlM')  TW  Paw'^p 

)  Fi  r  at  «T  ai  (  ?l  1  : 

nri.  rr,MS/bn:>  f^tpyi 

FI  OAT  PTA'(?n«  float  uTA|(71>,  float  WIM  (?])• 

F|  OAT  RT'''(?1)«  float  F(  OAT  (?1) 

)  t 

OOL  I  TFMP  F]  XFO 
nCl  <  PA  (fr('  )  « 
o  ( Sf! ,  )  • 

otFF  (F.OtSO)  • 

X (SO)  ♦ Y (505  , 

1 1  I  5  0  »  5  0  )  » 
iip(5n,5ni  « 

A  7  A  ( *=  0  I  » 

a  V  T  O  (  T  0  5  . 

APP  (30)1  FI.  OAT  t'TM! 

iiri  roivP  FMTP''  (  (*)i^LOST  =^T'''»  (•)^L  ’AT  PTN»FTXhn  pTK  * 

K  I  0  A  T  R  T  M  )  ; 
nn  (  T1  TNird)  5 
FTXrn  RTMJ 

nn  (  TFi  AO.  viflao,  I^laoi,  mflaoI)  '< »  T  ( 1 )  fxtcPMAl  : 

OCI  (  PCOPl«  POOP)  float 
nri.  OOF  Fi  OAT  BTO'J 

net  (  NO'-')  ftxFO  ptn? 
ori.  (  L  A  5*^  • ^  5P  » A  • A^  T  )  FfxFn 
noL  (  MST?f  VFTA)  FTX^O  PT^! 

ori.  (  5Tr7(4?),  PPX{47), 

«;TCA(4>'),  pSAIA?))  fiOaT  RI'iJ 
OM  FMOFTI  F  (FY5T‘‘'»  RFoTm; 

D|iT  pifTO(A)  L.TST  (IFM''  OF  OaTA  CApOS  -HAI  T  •  )  » 

OTOPt 
FKjn ; 

/#  OFT  A5P-B  OATA  A  'O  POPfTTO^I  * 

fi,r^  |T^T  (  MAW»(  AVT'MTM*  APP(Tl)  O'  Tlrl  To  N,AW))r'’®YI 
OFT  I  TST  f  AFIl.  F,  aF  -fi  )'~npY: 


r  A )  I 


Pege  Twent '/-Six-A 


i"n  >D  1 1  FW 


liP'-ir  |■■to  T  T  OMt;  ( M ft  I  rj )  t 


'•CT  I  r<;T  (  'JO«  MftO)r')PYI 
'■•^T  I  TST  f  ofo)  r^ovi 

fif”!  (  47_r'Ft.l»  A7oFi  l^XKPn  Pi  n?T  ►<!  J! 
n^'i  (  TftPoo,  rft<;r?A7)  Pin*T  ht'-i  (?])t 

neT  i.t«;t  f  A/OF|  1.  oroPl) 

pPT  I  r‘>r  (  xxa\  )  CPPYI 
r;i.- r  I  t«;t  r  '/\ipr<;p)  r^pvt 

P-PT  I  t«;t  (  ''5T?.  (  pp?(T)t  PTCP(T)  dm  isi  t")  nc  rp )  )  roPY ; 
pFT  I  TST  (  MPT4.  f  p^4(T),  STr4(T>  >  T  =  1  TO  MST4))CnDYI 

Gft  nft  t  a  (  c-iA  Ml  )  ; 

GFMSp:  r-PT  I.  IGT  (  '^APTf  MA*  TPl  A'il.  mFIAG)  COPY! 

"iWsM  A  W  t 
-VTTsAVT'^t 
wOPsAPR I 
O AO Afi  =  *  1  •  P ; 

rAl  I  «/SSrA'l(  AFTiP.  AP^'PL»  ■■’AST*  mA*  MO*  MAP.  ')*  IIP*  LAPP. 

Y  A  c  I  >  t 

TP  (  AS^>  Sfi  I  '^'ASR>  PO  TMFA’  l)0; 

PUT  FT^r  (•  ASP  ARRAY  STORA'^F  FYCFPPPO'T  (X(ln)*AJI 
stop; 

PUPS 

/«  r.r'ifRATF  a/tmiitmp  amo  oamGP  array  */ 

Piif  s  A7?  ♦  n,5  *  A7_f'r,l.  1  *  TMAR-l)  ♦  A7PPLU 
no  TsT  TO  I  ASRI  ” 
a7A ( T )  a  SUMI 

<;UMS  Pi|«  •»  (  '■■AR  -1>  A7__nPl  1  J 

PMP  ! 

PIIM  =  RAMPF  ♦  ^,S  *  PAMGF_nP|.  *  (MO-1  '  ♦  PCOPT  « 

no  T  =  1  TO  '’APR;  ” 

RA  f  T  )  =  cim: 

PI)'-<  =  SUA'  *  (MO-Tl  *PAMOF_'AF|  J 
FA'P ; 

/*  cibtniT  out  ASR-'I  nnwFP  dppfTVPP  OR  »FFlF'"TT  VTTY  */ 
on  T=  1  TO  I  ASRJ 
nri  u  =  T  TO  MASP! 

▼P  IMT.J)  >  A\/!n('IAW)  TMPM  oo: 

■KT.vll  = 

GO  TO  MFXTPl 
PMOT 

iKT.I)  s  PT  PI  (U  <  T  *>l)  .  AVTO,  maw,  APR)  I 
TFmPA  =  PTFl  (RA  (  J)  *oS4*'ST4.«;Tr4)  J 
TFmpP  =  .|(T,J)  ♦  TFmPaJ 

TP'APP  =  TFMPR  /TO.Oj 
TPUPO  -  in,  »«  (  TP'«P‘-')| 

U(T«J)  =(PA(J)**?*  tpmpp/xk-pu  *1  ,p  1  h*  <,>*1*1  ,Ffti 
'i(T*j)  =  IP,  «  1. ooTn(i'(I*n)« 

TP  I'^T.J)  O’f'.  T‘(PM  ufT«l)  =  11.; 

MFXTi-'J  PmO; 


f.n : 


OMP  J 


Pa^e  Twen  y-Seven 


:  iJi'or  liar  T  OMS  ( t '\i )  : 


I  r  t  A<;o : 

ir  I.ACP  >  30  thfm  I  =  30» 

PIlT  iTt;!  M  A<;p-a  '^ATa  f;} 

OUT  ';t<TP  rnjT  (  ♦  -p-».  (a/a(’>  n^'  t  =  l  to  i  aco)  ^ 

(Xr?),  A.  3n  F(&.0),  ^kTP)5 

on  jsi  TO  •jfAFki 

P((T  FtCfC  FOTT  (  PAfJ).  (  (ur».»)  no  fR]  TO  LA«=?)) 
30  F(4.0>«  FkTD)  : 

FvO ; 

LACO  =  LTFmp; 

no  j  =  1  TO  MAS»: 
no  T  r  1  TO  14*^0* 
n  (  T  .  J )  =  "  (  T  ,  j )  ; 

'l(I..J»  =  -1-»0I 

.IP(T«,I)  =-l?0i 
FNOI 
PND  J 

PUT  PAGE  I  1ST  (  •  •  )  t 

FL  AO  c  « 1 ‘G  t 

GFT  I  T?T  (  M£(  El  5 

nCL  (  AETiF*  ftENFl  )  FTXFn  BIN  (  il)« 
on  (  RANRE1/S0)« 

AZl (SO ) ) Float  qTNt 

ncL  wahaw  rtt  (1)  fvtkpma.i  tntt 

fiET  I.  TST  (  NP.  NA )  I 
OGL  (  W A (NA  »  WP1  , 

WAP  (  NA»  NIP)  , 

FTO  (  NA  «  NP)  . 

VAP(  NA»  NPj)  roMTPOLLEn  float  PTNI 
A(  LnCATF  WA  (NAf  MR)  » 

WAP (  NA«  NP) « 

STG  (  NA »  MP )  . 

VAP (  MA«  MP1 ( 

ncL  (  wvTn(3n)*  wppoO))  float  Rtn  external* 
nci  (  PANOiF^DFL  tntt  (0,07)  » 

AZ^riFL  TMIT  (0,07), 

A7REL  ♦ 

p» 

FUN)  Float  rtni 
nCL  Nw  FIXED  BIN  FXTFRNAL* 

OGL  (  WFTlE  ,  WFnFL)  ETXFO  PIM(31)J 
OGL  (  M‘;TART  ,  M,  NAl,  NA7,  NA7EN*  NRAN)  FIXFO 
OGL  WSSCAM  ENTRY  t  FIXED  «TN  (31).  FTXFO  RTN(31), 
FIXED  HTM,  FTXEO  MIN,  PTXPD  HTN,  F I XFD  PIN, 

(•,«)  F|  oat  RTN,  (*,#)  F|,OAT  PTN, 

FTXFD  PIN,  FTXFD  PTM) » 

NA3l  )  t 
ETO*n, n » 
v/APsn,('  I 


rnMPTLE^^ 


wCai  I  :  poor  f^PTro^S  (ma  i^) 
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waPsri ,  n : 

r;rT  |  wc-.ji-i  )  -ooyj 

r,FT  !  C'tfl  1  .  ►I'vri  C')py: 

Rft  I  i<;t  (  '^/_f^F(.,  /•■'fl  .prop)  tomyi 

OPT  I  T ST  ( xto )  ropy j 
OPT  I  TST  I  V^JOTSF)  < 

oft  l  rsr  (  mw.(  w\/Tn(T3)«  ppp{t3)  do  to  MW))rnpY« 

n  =1 : 

ONF :  on  TVsi  rn  MFi.Ft 

FXTRA:  oft  LTST  {  MST4PT.  *A,  TPl.ftO,  i-lFl,  AG )  HOPY ; 

OM  FMOFT)  F  (  SYSTM)  HFOINJJ 

OUT  SKT0(4)  I  TST  <  «PNn  data  C APOS-HALT ♦ ) » 

sTnpj 

FMDJ 

PaDAPs*0‘P« 

FALL  l«(SSCAAJ(i«i(FTLF,WFwFL,M<:TflPT.M.MAl  .  MA?,  WAP,  VaP, 

MAZFN,  NPAM)  I 

TF  VAZEN  >  so  (  URAM  >  SO  TMKN  DOJ 

OUT  EDIT  MWSW  APPAY  PTORAGE  KXCFFDEDM  fX(10),A)t 
STOP* 

Euni 

nCL  (  RAMGE  ,  A72)  FLOAT  qtM  EXTFRNAH 
TF  flag  TNfM  0O< 

C((Ms  a72  ♦  O.S  *A7«nEl  *  (MA7-1 )  ♦  AZ«FL» 
no  !al  TO  MAZFMJ 
a71  (  T  )  e  SIIMJ 

«:(iM  s  PllM  «  A7_OFl  *r-AZ->l)l 
FmO  { 

<;(IM  -  OflMOF  ♦  n.Pi  •  R '■.MGF__nFL  OfA'A]-))  *  OCDQt 

nn  1  =  1  TO  MO  a  A  ; 

PANOFl  (1)  =  FIJ'*! 

P((M  a  SUM  ♦  RAVJOF^OPL  *(NtM  -1)1 
FND  t  ” 

FA’O  ! 

FI  AG  a  tfUBj 

on.  OIFl  rrNTRy  (  float  RTN,  (*)  float  RTm,  FTXFO  rtm, 

(*)  FIOaT  BIN)  RFTJRMS  t  FLOAT  PIM)» 
on.  (  XKP)  F^OAT  UTmT 

nn.  (  XTM)  Float  otm; 

XTMsTMt 

no  Jal  TO  URAm; 
no  T=1  TP  MA/pMi 
p*  ramgf  1  ( .))  j 

TF  WADfT,.f)  >  .vVTO(MW)  THFf  on; 

,,AP  (  T  ,  J)  =OQ,  • 

GO  TO  NF.XTLI 

FMp; 

WAPIT,J)  =  PTFl  (WAP(T,J),  UVTO.NWfWPP) I 

tfmpa  e  wap(T,.i)  /in.n; 


CO-^PTL^-'R 


WCAli:  PPor  OPTIONiS  (MAIN)  t 
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WaP(1»J)  s  10,  *#(TrMPA)l 

wfiP  (ifj)  «  (rP**?  •  »(AP(T..n>  /y*fP)  •  i,fi«  •  3.a3  * 
WAP(TtJ)  *  10,  •LORI  0 (WAP ( T .J) ) J 
TP  WAP(I.J)  <  ?0  TMPM  WAP(T.J>  s  11. « 

VAP(T«J)  a  PTP1(  VAP(t*Jt«  WVTH,  KIM*  WPRM 
TCMP4  s  VAP(I.J)  /lO.Ot 
VaP(T.J)  a  1 0,*« (TFMPA) I 

VAPrT*j)  a  <  (o**?*VAP ( T •  jn /xwo) *1  .t.  1  H«:<, a3*  1  ,cf. t 
VAP  (T*J)a  10,  *  I  RRl 0 (VAP ( T »J)  '  ' 

FKin  I 


PKint 

Kl  a  NA7FMI 
KP  a  MRAMt 
DO  Jxl  TO  tfpt 

TFMPR  a  RAWGEI  ( J)  I 


/*  W<;p  RAWfiF  •/ 


DO  T 


001 


1  TO  W1  I 

TFMPA  a  A71(T)I  /*  WSP  A7IMt)TH  •/ 

0A|  I.  rOMVAP?  (RO.Tw*  TFMPR, TFMPA.  ARS*AAS)  I 
00  TCOM  a  1  To  mafRJ 
no  TON  a  1  TO  I  ASRf 

TF  AAS  <a  A7A(|A*;P)  A  AA^  >  A2A(1>  tWFKI 
7F  AAS<AAS-A7A(TONn  <»  A7FP  TmEM  DOI 
JAafOMI 

TF  AR^fAPR  -  RA(TrON>)  <a  OFP  THEN 
JP  a  TCONI 

MA(JA«.iP)  a  MAP(T*.IM 

VAP(T*.|)  I 
-  WA  {  JA..IR)  « 

WA  (.IA..IR)  <?0,  THFki  nTTF(j4,j») 

FMOl 

r.in  X 


not 


SI  I  A.  JO)  a 

0TFF(JA*J0)  a  0(JA.,IR) 
TF  n(jA«jP)<?r  ! 


T  1 


FnDi 

FMI' J 

FKini 

fno« 

FmP  t 

I  TFMP  a  I.  ASP* 

IF  LASP  >  30 
PilT  PA(;F  I  TST 
P()T  SKIP  FOTT 
LIST  (  • 
c;kTP  FOTT  (  • 


P|iT 

P|iT 


TMFN  I  ASR  a  TO* 

( *  "SRS^  data  *  )  I 

(t  mFom  VAIJIFS  •)  rX(20)»A)l 

A7THUTH  >..>« )  I 


«  (A7A(T) 

30  FC4,0) 


1 


1  TO 
FOTT 
TO 


nr  ja 
PIJT  SKIP 

no  I 
(X  (?)  .  F (A.  1  )  , 
FMOl 

oiiT  SFTP  FOTT  (* 


.  P-  f 

(  X  (?)  .A* 

HASP  « 

(  PA  (,  )  .  (  WA  (  T  .  J) 

ASP)  ) 

30  F (4,0) ,  SKIP) * 


00  T*1  TO 

.  SKIP)t 


I  ASH) ) 


RMc,  \/A|IIFS  ♦)  (  S*fIP  tXcPO).  A  »SKTP)» 


I 
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i*  rAli  t  »0'*»r  opt  T  o>'«;  (  «^a  |  m  )  j 


P.|T  «;t<TP  FOTT  (  •  -Q-  *,  (A/a(T)  :\n  Tsl  rr,  |  a  <;  "j  i  ) 

rxf?5.  A,  F(a,A'.  s<’P); 
nr.  .1=1  TO  MflSD ; 

r.i)T  SKIP  r.rTT  (  OAf.i).  (  STO(T..J)  i»0  1=1  to  i  a'I^M 
(X(P).  ►(6.1),  rn  F(4,0»,  Ck-TP); 

pmO  : 

X  T  vi=  T  M  : 

1)0  T=  1  TO  I 

n('  ,1=1  TO  ■uiflcpi 

TP  XTN*  =1  ThPv  TFmOa  =  10.  **  (II  (  T  .  J)  /  1  0  ,  )  t 

n.SF  TFmOa  =  1  0  ,**  (II  (  T  ,.J)  /  1  I)  ,  )  *  (  X  TM-n  J 

TPs»pp  =  WA(I.J)  /lO.O; 

TpMPf,  r  rpMPA  «  1  0  ,  *«  (  TF'*PR  )  : 

TfMCA  =  TFMOA  /  XTM» 

IF  TFMPA  <  =  n.O  TuP^i  TFmoa  =  1,aj 
U(T..J)  =  TFMPAt 

II  (  T  ,.1)  =  1  0  ,  '•>1,  roi  A  ( |.  (  T  ,  .J)  )  t 

Tp  yTn=i  thFni  tf*4i>a=  1  0  .  *<►  (IIP  (  r  ,.) )  /  1  n  , )  j 

flrf  tf.ipa  =  1 (MP  (  T .  J)  /  1  0,  )*(  X  T'4-.i  )  ; 

TFMop  =  sir,  (T.J)  /10.0S 

TFMpA  s  TFMPA  ♦  lA.n  (TFMPP): 

TFMPA  =  TFMPA  /  XJMt 

IIP  (T.J)  *  TFmbaI 

'ip(T,.i)  =  I'l.  »i  on  n  (IIP  ( T .,)) )  5 
e  'in ! 

F.Mn : 

PUT  <"KTP  fott  (  •  avFRagf  of  <;CAM<;t )  (,<  (?n)  .A)  t 

niiT  <;kTP  fOTT  (  »  -p-  »,  (a7A(T>  00  1  =  1  TO  LARR'> 

(X  (?)  ,  A,  30  F  14.0)  ,  Si'TP)  ! 
no  j=i  TO  maro! 

OUT  RKTP  FOTT  (  PA  (.1)  ♦  Ci  (  I  ,  I)  no  T  =  l  Tn  I«SP)) 

(X(?),  F(6,l),  10  F(4,n)»  RkTP)* 

FMO* 

PIIT  SKTP  FOTT  (  t  RMS  OF  SCAMS  •)  (X(?n),  A)« 

OIIT  FkIP  FOTT  (  •  -p-  •.  (A7A(T>  no  T  =  1  TO  lASP)) 

(X(?),  A.  3  0  F(4.0),  Si^TP)  I 
on  Jal  TO  MASPJ 

PUT  SKTP  FnTT(RA(.)),  (  iIB(T..M  on  1  =  1  TO  |_aSP)) 

(X(?),  F(6,l).  10  F(4,n),  SKTP) « 

fmOJ 

piiT  PAPF  FOTT  (*  FA4  PAUAM  MTMUS  WSP'  ) 

(  X  (1  0)  ,  A  )  » 

no  .1=  1  TO  -lasRi 

PilT  FkTP  FOTT  (  PA  ( J)  •  (OTCF  (  I  ,  J)  00  T  =  1  TO  LASP)) 

( V ( ?)  ,F (6, 1  )  ,  in  F(4,0),  SKTP)! 

fmO  t 

pi|T  OARF  FOIT  (  *  '”0RPF|.  A  TTOM  TM  A'TMlITH  *)  (X(ln),  v)l 

nr  .Jsl  TO  MASPI 
no  T=1  TO  I  asp: 


T,_c 


p  =  or  (apt|C<*,  n 
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•«;TAPT.*-'»M.'M»  '1*7,  WAS*  Vfl>S. 


I  Tf; TTv'^ 


po'^r  FMfiir.  apt  » « •^|  a  i  .  ma7,  vap*;# 

M  A  /  r  M «  M  p  A  M )  ; 

nci  t  vAof;(o,«))  PioAT  hiaii 

l-x-l  t  WVTHfTDt  Kinaf3li))  C|  P  T  A‘  ^XlpP-iAL  S 

(  AFTlF  )  FTKPn  PfNfin: 
nri_  t  ox(30)*  PYCiOn  FiOAT  ptm* 
n^i  Mr>  p  T  XF  )  m  V  : 

nrt  papam  RTT(i)  fxtF»mai,  tvTT 
nri  A  xi  ftyFT  fxtfhmai  : 

nr(_  (  WA(?A,  s«))  F|OAT  ptm; 
nr;,  t  vaP(?6*  =;b)^  fL'^at  pT'j; 
nrL  (  VAl,  MA/)  e-TXPn  HTM| 
nrL  (  a7F'>iTm  tmt  T  f  0  i  , 

pfmtm  TMTr(O))  Frxrn  PT-t; 

nri.  f F()m,cu'-<i  )  Ft^oAT  MTv: 

nr;,  ctfi  fmtsv  (Fi  paT  ptm,  <*)  fioav  ni^i,  FTxm  btm. 

(•)  float  PIM)  PFTilRAP  fFI  PAT  MTM)  » 
nrL  rHAML  r-tAR<n  ExTfrmm  t 
nFL  tfmoi  Float  nfMt 

DC(  r-jpc.  ,MnMTH«r)AY.  YFAO»woM*p(?F,wT_Typt  .'^CAm^MO.Pt  )  FTXFO  utm: 

nri.  C^)  ftxf")  him  ” 

OOI  VI  (  154?)  FLOAT  PT’'-; 
nri  VTOFO  RTT(1P0001  \'ARt 
^C|  T,*.PF  FTLF  tmo'IT  » 
on.  f  fnfti  F)  FTXt^n 

nn  (  RA‘IRF  .  A7?)  FI  OAT  PTM  FXTfPNAI  5 

on  TapPM  FT|  F  tmoiitj 
ori  (  AFTIF*'  FjxFf)  PiMfin; 

nri  •'PCFvp  F  VTPY  (F  T|  c- .  FTVFr)  P  T  Ni  (  3  1  )  ,  CNAP(]7m: 

nri  rrnF  pyx-n  SF^'SF„o3n  chap  (1P)J 

nri  (MCiiNi.X'MFC,  TCnnF)  ftxfo  nTM(31)J 
nri.  ( T 1 )  r-TXFn  py^'  FXTFPrjai  i 
nri  f  K17.  <)fl.  yx.  T *  I ,  yx7. 

T7,  yp,  J,  H,  M-lt  KlP,K?,  K?n,  <3,  K4*  XF.  «7,  I  IP, I  4, 

I  H,  ypAM,  tSTAPT*  I'fPflAVPt  yAV3.  TFX  lyMyTlAl 

F  T  YFr  P  T  NJ  X 

Of;  OPFA'Fyi  FMTOY  (C^-IAO  (F)  ,Fy  xf:)  RyAj  (  3 1  )  »  f  y  yFO  PTMiyin* 

'''"1.  (  A7» 

Kio_PT^,  /•  MO  nr  vynfo  oys  pojmtfo  fxiT  */ 

“O  woROF,  /»  VO  OF  coMt^ijIFw  t.inpiip  nm  «/ 


kio_pt^  « 

«  o__woRnF . 

y  o  p  n , 
••r «; a n  . 

•  C  l.i  F  F  3  , 

/"oFl  , 

;;  a»  i-c  _nFi  • 

A  X  MF  AM  t 


/•  MFMOJY  STAPy  AnnPF5<; 

/»  ^'o  nF  i;'4(Fi-ps  In  a  wyvOTw 
/*  a  7  TMi  iTH  nr  I  y  A 
/*  RaMr.F  nri  ta  »/ 

/*  .*ay  vxiTMnRY  ftO'iRFFp 


r 


ro  IP  fL»- 


Pace  T1  irtv-Threp-A 

wSsTAm:  (Arri>-rt.  K  ^  ►  I  I  '•  f  "Q  T  Af  T  ,  M  ,  I,  ^  1  ,  IVJA7,  MA‘n.  \J  M 


“'^Lf'^SWFrP  /*  MTi  e  of  j  P 

)  e-T  xFO  ( *^1  ’  : 

nri  f  •^ftarT)  Fixm  ^»*i; 

/«  mctapt  Ttjc  Rcrrtpn  •/ 

/*  M  T«;  TPF  FTTP  RPrrlpr^  «/ 
vo_PTFs  J 

Di(su(\,Tr>: 

PV  =  WPP  » 

^|PS^.W^ 

ON  rMPitTif-  f<;v«;TN^  PFfif  t 

P(lT  ^KTO(<*i  -TST  (»CM!'  OATA  c  AP’''*;-HAt.  T  »  )  t 

FTOPI 

r‘'OJ 

ON  tOAN«;utt  (T4JF)  qc»;t^4; 

CAij  TOFFK<;(TAPF«CnnF.c*  NScwnp.>)  : 
jf  ronF  =  1  irnoFslI  rnoFsA  tufm  ';o  to  i.  <; 

PiiT  FtfTPMl  I.  T  (  ♦  TPAM^'v^yT  FPp  rOOF*  •  H  CO'^F  I  1  '  • 

•SrNFF.  T*?  •  I  I I 

ttM  TO  I  0 1 
fnd; 

o*'  TOAMftM|T  (TapF(^)  MPOTNl 

TALL  TPSFm‘;(  TAPFV*  r'^OF.  FFN^FWOPO)  : 

TF  ropFsi  (ronF*3i  roncs4  twpm  fio  to  l3* 

pllT  SKTPfl)  I  I  ST  (  •  TOAm^v-jx  coo  roOFe*  I  ICODCM  •MPCCfcMO'lMs'  ♦ 

*^FM^F'-OPn  TS  •  I  I  ^FNSF  »«I0PI>)  I 

oo  TO  I  OJ 
pwn: 

OM  Fr'OCTI.F  (TAOFwi)q^GT''': 

PIIT  FK(0  (?)  LTST  (tcMO  OF  F  TLt  '••FArwEnM: 

STOP  » 
p»>p  t 

ON  FNOCTLF  (TAPC)  hfotv; 

PIIT  OKI®  (?)  I  tST  (icma  nr  riLfc  ‘^FACmEDM: 

STOP  t 

FNOi 

ACTl.C  »  AFTLF»** 

•f;:  TF  AFTL.F  >  fmFILF  than 
no  t 

PMT  RtcTP  (?)  I.  TST  (  *  All  DFrOPOq  WFAO  *)« 

pFtiirm; 

fmoi 

TF  PAOAP  THFm  001 
ri  OSF  KTI.  F  (TAPE)  S 

Tail  OPEnctl*  *TAPr»,  ACTir,  tCoOF)! 

TF  TOfiOF  =1  TmEN 

oo* 

PIIT  sktpi'-')  list  mmaO  opaO  Im  opfm  cTlf  -HaiT*)« 

OTOP  t 
cmoi 
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p  r f- 


ODor  (aftipk*  f-vf-it!-,  ■<«;  T  ap  t  ,  ,  ki  a  i  ,  ma7.  '«/as» 


TC-  Th?m  Tl=lt 

T nwr  =n  ; 

T'^i'"=l  ; 

TO  (IMF; 

F^ln  t 

tlocf  FTI  F  (  Taofw  ^  I 

''All  ."1PCmF^_(  iTAPF.!'.  AP'TIF.  TFri^t'S 
tf  TFnnF  =1  THF^! 

''.o ; 

OUT  FKTDi?)  ltst  (ti'An  'jfao  T^(  npf'i  ftif  -malt*); 
<;tod  ! 

r  \|  D  ] 

TF  TFL  Art  THF  N/  T  1  =  1  t 
T  ooFsO  ; 

T  «;  A  M  r  ■(  ! 

nn  v^oTSF  f^oar  mtm  •xtf^njai.: 

n''F;  TF  OaHA-?  TmFM  OFaD  F  Tl  F  (  T  APF  )  1  \J  '  O  n/THFA): 

FL^F  pfaO  FTl  f  (  TA3FU)  ImT'1  (  v'Tofo): 

TF  Ti<  mftadt  tmfm  f;n  r-i  )_?; 
vn]?:  nn  tfI  to  A'-\_dtsj 
Kl  s  ?n«(T-1)  ♦  I: 

K<J  =  ?r'*(T-1)  ♦  ns 

TF  pahap  *  c-<ANL=*n  tmfm  \/i(t)=  s'ih's  tb  ( V  t  •  <  n  n )  s 

FinP"  VUT)  =  S"H'?TP  (  V  T  )FO»  K?,  }'^)t 

Tf-'pi  =  VT  (  T  )  ; 

TEMPI  -  TFMPl  -  V'lOTFFt 
\;  M  T  )  =  tfmp  1  S 

frio ! 

aZ=  FI iPFTB  ( V Tf>FT) »??<■),  1  ? )  «•  ,nA79i?? 

FiiRFTP  (V  nFO,?!  T  ,  1  )  ; 

Ain_pnOnF  =  Fl!PFTJ(VTnFn,lF7»P)  I 
tf  tdpf  -n  Them  not 

TF  aUih^s)  '40  TO  (  ; 

on  TO  T  OF  0  ; 

F»'r, ; 

TF  Tpoon  =  '■j')_/jnpnc;  fL  'IIImsT  TmFm  RO  T<‘  TOFO; 

F|  SF  on  TO  I  ; 

l(,c.n>  ..rOAO  =  FiiHFTPT  VfOFO,  ?)■’♦  h) 

T  F  AT  r  It 

RAf  r-F  =  Fll'JFTB  (  V  T  TFo  ,  ^  { 

Hk,Fr  CmO«;to  (  VT  nFO,  )  .R)  t 
''0^'TH  =  F(  iPFTP  (  V  T  OFo  ,  T  n  .  o )  : 

P7V=  ‘.Ht-FTP  f  V  T  :)F0 . 1  <- i  I 
Yt-AOs  e  I'BF  TO  (  V  T  DEO  ,  ?  T  ,  B  )  ; 

RP"  =  E|iPFTP(VTpFn,‘,7,A,); 

PwF  =  fi'uctr  ( v T nEo « 7:1 .  n  )  s 

V1_*YPF  =  FHUR  TO  (  V  T  oFo  f  ♦  *< )  : 

Fr/'M_',|o  r  CHRF  TO  (  V  Torn  .  !  Ffr  .  1  )  S 

ST  =  F:  imftm  r  V' T  "'Fn ,  1  K?  , />  >  ; 


V  A  PS 


roMp y 


vCQf^M•  non-'  (iPTIfM 


1  ).!- 
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••<;T  ftOT  ,  M  ,  Nft  1  ,  \)a7.  wa«;. 


Tunon  =  Mn_i»/npne;t 

■^'AVtiTAn  I-  «.|ipQTP  f  wT  'FO  .  IPP,  1  ; 

A.o^cmrcD  -  1=,?7  /  (mO_wOO<a<;  ny/-?',-. 

npL  Pi^AAT 

=  o.a; 

TFXsr  ; 

on.  (  Tt  CTT  IMTT  (  f'  ^  . 

MCI  AR  TMITtni)  FT'/FO  -IMt 
MCI  AP  =  0  t 
Ti.FPTrO  t 
T  AV^=n  t 

7c  npM  =n,n  thpm  oo: 

OUT  CKTp(p)  LTST  If  poM=n»>t 
OPM  =  ^,0t 

fMO  t 

ncr,__i^wr(rp  =  ■^AO,  /  oP»  /  lo'^n,: 

A?^^  =  ncfi_",WFFP  • 

«;  T  =  <  )  I  ** 

PAMRr_nPL  =  1,  /STt 

MTt?F__«:wFro  =  Mn_.  opo<;  •  oangf^dFL: 

oaMPF^IsTOp  =  panRF  ♦"’M|LF_SWFrp  t 
A7^«;T0P  s  a/  ♦  A^OFI  I  ” 
a7?  *  A7t 

Pi'T  <;kTP  FPTTCPErooh  Nr,*«T1)  fA»FM0*0))J 

PUT  <;kTP<?)  FDT  T  <  *  TTmf  :  *  »MAMTM  J  »  »M0NTM, ‘nAY  :  •  •OAY  ♦ 

«  yfap:  •  ♦vfad.  *oPw: »  «kp***  'pof  i »  .opf» 

♦  viopo  TYt’F  j »  ,wt_typf,  'scftMS  •  f<;rAM^MO, 

•SA'^Pl  T'ifi  TMTCPVAl  •  .SI  .  *0AMPF:  »  ,7aMRF. 

•A'O,  Of  rOWOUTFP  WOOOc/SwFPP  •  *Mn_'*'nPi)S  ♦  '  A?  TMIJTM  !  »  •  A7  ) 

(1(X(?)  .A.K(in,?)  »xT?)  .A.FdO.P)  .X(?>  .A.PdO,?^  , 
X(?).A,F(ln,PMM 
PI)T  PARC  COTT  ( 

I  MfKOPY  CTA^T  AOOPCSS  s  •»  mc<;ao, 

I  Mp  OF  St>'PP3  =  MO^SWFFO. 

•  OFROFF  PXTFviT  *.  0Cr”*»WPCO, 

•  RTOP  PA^'RF  PaMof^STOO, 

t  mtlfc;  ofo  swPfPt.  mti_f__SwFPo, 

•  OAA'RP  _  OFI  »,  OAMOC_OF|  , 

t  A7  STOP  »,  A7_STnp)  {  X(?0).  A*  F(S)»SKTP.  X(?0>, 

A,  PfS.o),  skt^T  y(P0).  A,  P(s.7),  SKIP.  y(?0),  a,  F(S«0), 

PK  TP,  X  (  ?fM  ,  A  .  F  (R,?)  , 

RkTp,  t (?n)  ,  A. 
t  DOF  s T : 

TST aRTcQ  s 
\/ap  =  n,n; 
wAsn , 0 1 

/\-»FfjTH  r  mO^SYPPP  /ma?« 

I  fmtm  s  ywoD'i  *T/pj 

rJF^iyH  =  (CmT-H  /MAIt 
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'dtl^ 


L 


'./SSrflM;  Dr,0':  (^^T(r<.,  FMf-Jl.i-,  m  ,  m  a  )  .  ^l  ft  /  , 


T  ft  V?=  I  r*  tm ; 

7  K‘^=)  ! 

yp  MC|Af?  =  0  T^F^J 

T  A  \  /  ■?  =  ■ )  l 

71  c-FT  =  -  '\JO_'^WK  Km  *1  p\TKi  _7STA«TJ 

r>  -^0  o  ?  n  ; 

KAiO : 

TA\/.^  =  TL^FT-  4.  TKyJ 

Tl  FKT  =  l';'7/S  ♦  71.  f^FT  -  MO_t;  jfro  *iKMTM| 

MKI.ftC  =0S 

!•<'=•  =1  ; 

7F  TFX  -.  =  ri  t-(knj  nr>| 

TK  TFy  ■>  I  FixiTM  T'-'fN  n<'j 

P‘IT  1  T^'T  (  TFX  )  : 

TFX  =')! 
r;n  Tn  upg* 

FMn ! 

nn  k7=l  Tr  r-v{ 

TFMD  1  -  '/ 1  { ><  7  )  ; 

quw?  =  SI  IMP  *  tfmo  1  t 

SIIM1  s  SIIMl  *.  TFMPl  »«p; 

FmH  5 

[F  Tx  >sn  t«fm  no* 

oi)T  CK  TO  FOT  T  (  IX.  »  »  ’  ) 

(X  MO)  ♦  F(iS«0)  »ft)  : 
r;n  to  opnj 
FNP  5 

V  A  {  1  ,  T  X  )  r  S  JM?  /A'M  J 
VAO(l,TX)  =  SliMlJ 
T‘<^  =  K7f 
FmO  5 

0?n;  nn  T7=1  to  Ain_SWFFOj 
X  v=  M  ♦TflV  T)  /'lAl  *1  ; 
on  T  =1  TO  LEMTh  by  NAlJ 

(T*TDV0)  >1  Fmth  thfm  on  TO 
L1«=T* 

(4  =  TtMAI  -1  ! 

SUM  =n,; 

=  n.ot 

no  I  s  =  ij  A  TO  14: 
kS  r  IkS  ♦  TOTftPTI 

TF  kK  TMKAf  on  ro  n?i  I 

Tf  (14*  TAyO)  >  I  FNTm  Tmfm  no  TO 
TFMP  1  =  '/  1  (  ‘'’F  )  t 

OHM  r  IM  *  TF  IP  1  t 

SUM)  s  SUM)  ♦  TFMoi  **p: 

T  k  =-  =  T  kS  ♦  1  : 

F  mIT  ; 

T  ^  I  y  >  s  n  T  i-i "  M  n  n  : 


WAS  . 


Ff-> 


: 

1  ! 


r)?6 : 

OPS 
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VSSrAN':  Opor  (AFTIf.#,  PMf-lIK,  MCTftOT^'<*^•Al  .  MA7,  W  A 'i  f  VAPS. 


PJT  SKTP  E  'TT  f  TH,  ) 

(vdO't 

TO  OSA,: 

Fmo  j 

|.'A  r  T7,  Jy)  =  5ll»«/‘tAl  I 

\/aP(T7,  Tv^  =  ^1)41  ; 

T  X  =  T  X  ♦  ]» 

!  pmh: 

TA''P=0j 

FMft 

oo  TXI  =1  TO  RENThi 

no  Jrl  TO  M0_SWFFo  py  AJA7t 

KP  =  j: 

kA  s  J*MA7  -II 
TX?  =  .J/MA7  *11 

siioso.o: 

SII**!  =  n,o« 

Do  =K3  TO  IT  4 1 

TF  Jl  >NO„S^PPP  TuPm  on  TO  0?S| 

SUM  s  sumTwa  (.11 ,  T X 1 )  I 

StlMl  a  SlIMl  ♦  VAO(  ,)1«  T<l)t 

fnDi 

WAHT?.  Txn  a  SIIM/MA/I 

SUM!  s  /(NA1*NA71  -  WA(TX?»  T*!)  •*P| 

TP  SiiMi  <0,0  THFN  noi 

VAP(  IXP,  TXT)  sO.OI 
Oo  TO  Q7SI 
PMO  ; 

\/AP(  TX?,  Txi)  r  SorTcsiimim 
T  X?  r  T  X?  *11 
!  PNO I 
F^'ni 

NA7FN  =  TS4V  *  A7PmTh: 

MPAM  =  PFMTHI 

K)7  a  (ISAM  -1)  <»A7FNTH  *11 

«'!«  =  <17  ♦  A7E»JTM  -i;  ■ 

k-1  o  =  1  I 

no  IRaKlT  TO  Kl«; 
no  1=1  TO  REMTMI 

I-'AS  (7R.  T)  a  waIkIM,  T)t 

PMO  I 

K-IP  =  Kl9  ♦!  I 

P  N' .  1  I 

ton  r  11 

Of'  TO  a  KIT  To  Kisi 

on  Tal  TO  pPvlTMl 

VAPS  (TP.  T)  a  WAOFKEn,  T)t 

F^ni 

KPOs  K?0  ♦ 1  I 
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)  ■v.p  T 


«/Ssr4»j;  ppnr  <  *fri  f  ji  ^  ,  •'^:T/*.  'jt.vi.nui,  m/*,  /  ,  wa'?* 


f'jO : 

r<^rfCT  =n: 

TP  TifFT  =n  thF^!  r;n  T-'  I 
TP  Ti.fpt  >  i_='MTH  THF*^  nn; 

no  T  =  '.  TO  _fmtt  rv  ;'Al  ! 
r  t  =  T/M/iT  *1: 

I  1 fi  =  T  : 

L4  =  T+M^l  -1» 

«;iiM=n .  0 » 

c"M  1  =  n ,  j 

on  L  <^r(  1  «  TO  L  4  ! 

K=:=7K'=«; 

TP  isi^  thF''!  on  T-'  pi'3; 

TP  I  I  FMT-t  ThF,-..  OO  to  OI't; 

TFK<P  1  =  v  1  <  <P  )  ! 

Piin  r  qi IP  *  r Fwm  1 


CM*'  1 

=  S  H  H 

1  ♦ 

TPAOI  *«?; 

T  Y  S  S 

TKP 

♦  1 : 

p '  1  n  } 

TP  T  y 

Twpm 

not 

OJT  PKTP  E'tTT  (  T'<« 

( y  (1  0 ) «  p  ( F. ,  0  T ,  a  >  I 
i?n  TO  PIP; 

Fvn ; 

i./a  ( 1  .  TXT  s  /  ma  1 ; 

VAP  ( 1  .  fy  )  =  'ii'MT  : 
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FAA  RADAR  GAIN  PATTERNS 


The  antenna  gain  patterns  for  the  ARSR-lD,  ARSR-2/1E, 
;!^SR-3,  ASR-8,  and  ASR-4/5/6/7  are  shown  in  Figures  B-1  through 
E-5.  The  curves  are  shown  plotted  as  gain  versus  relative  eleva 
t ion  angle.  The  actual  gain  at  a  given  elevation  angle  depends 
cn  the  positioning  of  the  antenna  in  elevation.  The  plotted 
c  irves  represent  data  points  read  at  .5“  elevation  ir.crements 
f  -ora  curves  furnished  by  the  FAA. 
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Figure  B 
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-1.  Antenna  gain  pattern  for  ARSR-lD  versus 
relative  elevation  angle. 
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Figure  B- 
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Elevation  angle  (degrees) 


Antenna  gain  pattern  for  ARSR-2/1E  versus 
relative  elevation  angle. 
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gure  B- 


Antenna  gain  pattern  for  ASR-8  versus  relative 
elevation  angle.  High  and  low  beams  shown. 


ATC  ENRO’JTE  WEATHER  DETECTION 


This  appendix  is  a  report  written  hy  W.  Goodchild  of 
K  \FEC  which  documents  the  weather  measurement  experiment  performed 
i  1  Oklahoma  during  September  and  October  1977, 
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DISCUSSION 


From  late  August  1977  to  November  1,  1977,  sinwltaneous  comparative  radar 
weather  data  was  collected  by  NAl'EC  and  National  Severe  Storm  Laboratory 
(NSSL)  personnel.  This  data  was  collected  from  a  typical  WSR-57  radar  which 
is  used  extensively  by  National  Weather  Service  (NWS)  for  severe  weather 
detection,  etc.,  and  from  a  typical  ARSR-1  radar  which  is  used  by  FAA  primarily 
for  aircraft  surveillance.  The  W5R-57  was  sitcated  at  NSSL  in  Norman, 

Oklahoma  and  the  ARSR-1  was  located  at  FAA  Aeronautical  Center  in 
Oklahoma  City.  The  Aeronautical  Center  is  located  approximate..y  12  nautical 
miles  northwest  of  NSSL. 

Both  the  WSR-57  and  ARSR-1  radars  were  operated  with  establish.'d  radar 
parameters.  A  dual  channel  video  recorder  was  added  to  the  instrumentation 
already  available  at  both  sites  for  simultaneous  data  collection.  Channel  1 
of  both  recorders  recorded  the  intensity  of  the  weather  (log  video  output  of 
the  radar) .  The  log  video  output  of  the  ARSR-1  was  recorded  at  the  output 
of  the  Cardion  modification.  The  log  video  output  of  both  radars  had  at  least 
65  dB  of  linear  dynamic  range.  Channel  2  of  the  recorder  at  NSSL  recorded 
the  elevation  of  the  WSR-57  antenna.  Tnis  information  was  recorded  so  that 
the  vertical  profile  of  weather  cells  could  be  determined.  On  channel  2 
of  the  recorder  at  the  ARSR-1  site,  the  moving  target  indicator  (im)  video 
was  recorded.  Tire  dynamic  range  of  the  MTI  video  was  25  dB.  The  auxiliary 
channels  of  both  recorders  recorded  time  (for  simultaneous  comparisons) , 
range,  and  azimuth  information.  Photographs  were  also  taken  at  both  sites 
as  secondary  data.  At  NSSL,  the  output  of  the  Video  Integrator  Processor 
(VIP)  was  photographed.  This  provided  six  levels  of  weather  intensity 
information  on  photographs.  At  the  ARSR-1  site,  photographs  were  also  taken 
but  since  no  VIP  was  available  there  was  no  level  information  on  the  photo¬ 
graphs.  Figures  1  and  2  are  block  diagram  of  data  collection  configurations. 

Approximately  20  hours  of  weather  data  were  gathered  at  both  sites.  Weather 
information  was  recorded  on  a  total  of  six  different  days.  The  weather  data 
that  was  collected  was  observed  to  be  from  different  types  of  precipitation 
Structures,  including  isolated  air  mass  showers,  squall  lines  accompanying 
a  front,  thunderstorms,  and  stratiform  rain. 

In  an  attempt  to  determine  if  the  ARSR-1  can  be  used  simultaneously  for 
weather  and  aircraft  detection  without  degradation  to  either,  the  ARSR-1 
was  operated  in  several  different  modes  during  the  weather  data  collection. 

The  weather  data  collected  under  these  different  conditions  was  recorded  for 
a  comparison  with  the  weather  data  collected  simultaneously  at  NSSL. 

Standard  calibration  procedures  were  performed  on  both  radar  systems  throughout 
the  teat  period.  Sunstrooe  data  -as  collected  to  permit  accurate  registration 
of  tlio  data  from  the  two  radar  sites.  Receiver  sensitivity  data  was  collected 
using  the  video  tape  recorders  to  provide  accurate  calibration  of  the  video 
radar  data  t.ipes. 
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Upon  complecicn  of  the  data  collection,  the  video  recordiitt’s  were  returned  to 
NAFEC  for  digitizing.  The  video  recordings  were  played  bacl-  into  a  dual 
channel  Video  Quantizer  and  i'^ecorder  (V(}K)  with  the  output  being  recorded 
on  a  high  speed  digital  recorder.  These  digital  recordings  were  then  :3enL 
to  the  Applied  Physics  Laboratory  (APL)  for  analysis.  The  analog  video 
recordings  were  also  used  at  NAFEC  to  determine  the  capability  of  the  Weather 
and  Fixed  Map  Unit  (WFMU)  to  accurately  depict  weather  under  the  phase  2 
portion  of  this  project. 

TESTS  CONDUCTED 

GENERAL 

Several  radar  parameters  greatly  effect  the  display  of  weather  versus  aircraft 
targets.  Antenna  polarization,  MTl  velocity  response,  sensitivity  time 
control,  and  weather  detection  over  ground  clutter  were  each  the  subject  of 
Individual  and  combined  data  collection  efforts  under  this  project.  Table  1 
lists  the  various  test  combinations  for  which  data  was  collected.  Each  of 
the  above  parameters/tests  is  discussed  below.  A  further  chronological 
breakdown  of  data  taken  and  available  for  analysis  is  given  in  table  2. 

Comparative  WSR-57  weather  data  was  collected  for  each  of  the  ARSR-1  data 
segments  listed  in  table  2  and  is  retrievable  using  the  same  data/tape/time 
information. 

ANTFMA  POLARIZATION 

This  test  was  performed  for  determination  of  the  effects  of  polarization  on 
weather  detection  and  to  investigate  use  of  the  orthogonal  circular  polar¬ 
ization  (OCP)  signal  to  provide  dedicated  weather  information. 

The  effects  of  antenna  polarization  on  aircraft  and  weather  detection  are 
well  known.  If  weather  conditions  are  present,  it  is  common  practice  to 
use  circular  polarization  (CP)  which  will  reduce  the  weather  signal  level 
by  cpproximacely  16  dB  (the  amount  of  decrease  depends  upon  the  sphericity 
of  the  rain  drops)  while  to  a  lesser  degree  (2  to  6  dB)  reducing  the  signal 
returns  from  aircraft.  The  level  and  real  extent  of  weather  returns  (desired 
for  weather  formatting  and  display)  thus  become  unknown  when  using  CP. 

However,  it  is  possible  in  the  case  of  the  ARSR-1  polarizer  to  use  the  OCP 
(re\erse  polarized)  signal  component  to  provide  dedicated  weather  information 
of  the  same  amplitude  as  that  available  using  linear  polarization. 

Since  CP  was  originally  derived  to  eliminate  weather  from  the  radar  display, 
the  OCP  signal  was  dissipated  into  f  „  .ched  termination.  Use  of  the  OCP 
signal  requires  that  a  low  power  RF  path  be  provided  to  transmit  the  OCP 
8ign.il  from  the  polarizer  to  the  radar  receiver.  This  path  Includes  an 
exCt.i  channel  in  the  rotary  joint. 


2 


Page  Six-C 


SENSITIVITY  TIME  CONTROL 

This  test  was  designed  to  determine  the  optimum  STC  characteristic  for  the 
detection  of  aircraft  versus  weather.  Two  differe  nt  STC  values  were 
used.  They  were  the  inverse  of  the  second  power  of  the  range  (R“2)  extending 
to  80  nautical  miles  (nmi)  and  the  Inverse  of  the  forth  power  of  the  range 
(R~^)  extending  to  50  nmi.  An  R~^  curve  provides  range  normalization  for 
discrete  targets  signals  since  their  signal  power  decreases  in  a  forth  power 
relationship  with  range.  However,  since  weather  exhibits  an  antenna  beam 
filling  characteristic,  its  return  signal  power  decreases  in  a  second  power 
manner  with  range.  Thus,  an  inverse  fourth  power  curve,  while  providing 
optimum  discrete  target  range  normalization  and  detection  '"crly  attenuates 
weather  signals.  A  second  power  curve  while  providing  optinrjm  wenther  detection 
doesn’t  correctly  normalize  aircraft  signals  and,  in  addition,  allows  excessive 
ground  clutter  signals  into  the  radar  receiver/processor.  ITiis  excessive 
input  of  spurious  signals  will  ultimately  result  in  a  higher  false  alarm  rate 
into  the  common  digitizer. 

The  R“^  maximum  range  value  (80  nmi)  was  chosen  due  to  the  effect  of  the 
curvature  of  the  earth  on  the  detection  of  radar  signals.  That  is,  while 
in  some  cases  the  most  intense  portion  of  a  weather  cell  is  close  to  the 
ground,  for  a  weather  cell  at  100  nmi  only  those  portions  above  6000  ft 
elevation  could  be  observed  with  the  radar.  This  problem  gets  worse  with 
increased  range  making  such  long  range  weather  data  unreliable.  Therefore, 
only  the  data  within  approximately  100  nmi  of  the  radar  site  was  used  and 
the  STC  curve  was  adjusted  accordingly.  The  R~^  curve  is  representative 
of  those  in  general  use  in  the  FAA. 


An  additional  feature  of  this  data  collection  sequence  is  that  the  R”2  curve 
for  weather  extends  further  in  range  than  the  typical  R~^  curve.  This  results 
in  a  crossover  phenomenon  resulting  in  the  R~2  curve  having  greater  attenuation 
tnan  the  R“^  curve  at  long  ranges  and  less  at  short  ranges.  Table  ]  lists 
the  various  tests  conducted  to  collect  data  relative  to  the  two  STC  curves 
and  other  pertinent  radar  parameters. 


SIGNAL  ATTENTUATION 


Tlie  next  category  of  data  listed  in  table  1  relates  to  signal  attentuation . 

This  test  was  designed  to  enchance  detection  of  weather  over  ground  clutter. 

MTl  systems  in  use  with  current  FAA  radars  utilize  limiting  of  the  intermediate 
fiequency  signals  prior  to  phase  detection.  This  is  done  to  keep  the  signals 
w  thin  the  djmamic  range  of  the  cancellers  thereby  keeping  ground  clutter 
rt  sldue  to  a  minimum.  However,  this  limiting  also  reduces  the  amount  of 
wi  .ither  available  for  display.  By  attenuating  the  radar  receiver  signals 
pi  lor  to  limiting,  adcitional  useful  weather  information  may  be  available. 

Til*  te.sts  listed  in  table  1  show  the  various  configurations  tested  and  the 
anount  of  attenuation  used  to  provide  the  desired  data  for  analysis. 
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MTI  VELOCITY  RESPONSE 

The  last  category  of  data  in  table  i  was  taken  with  respect  to  system  MTI 
velocity  response.  The  Maximum  setting  shown  provided  the  greatest  possible 
narrowing  of  the  MTI  fll-.er  low  radial  velocity  band  reject  notch.  This 
provided  the  maximum  amount  of  weather  information  for  analysis  while  still 
rejecting  ground  clutter.  Data  was  collected  for  analysis  for  the  related 
configurations  shown  in  lable  1. 

DATA  FORWARDED  TO  API.  FO T  ANALYSIS 

The  above  lists  of  data  .;ollected  were  forwarded  to  APL  along  with  the 
weather  cell  photographs  taken  at  the  time  of  data  collection.  Based  on  this 
Information,  the  weather  data  windows  shown  in  table  3  were  selected  for 
digitising  at  NAFEC.  The  data  reduction  configuration  shown  in  figure  3 
was  used  to  provide  the  desired  digital  tapes.  The  data  shown  in  table  3 
were  forwarded  to  APL.  Calibration  data  from  both  radars,  azimuth  correction 
(offset)  data,  and  ground  clutter  levels  with  STC  1  (R-4)  and  STC  2  (Pv-2) 
were  furnished. 

RESULTS 

1.  Seven  comparative  ARSR-l/WSR-57  weather  window  data  tape  pairs  were 
forwarded  to  APL  for  analysis  along  with  the  necessary  calibration  Information. 

2«  Video  data  tapes  were  made  for  studies  of  the  weather  quantization  and 
display  capabilities  of  the  WFMU  at  NAFEC. 
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This  evaluation  was  conducted  In  response  to  9550-AAT-100-33  "Provide  Improved 
Display  of  Weather  Data  on  Radar  Displays,"  dated  August  29,  1975;  Subprogram 
021-241,  sponsored  by  ARD-243,  Mr.  Kenneth  Coonley.  This  project  number  is 
021-241-130,  the  NAFEC  Program  Manager  is  Ronald  Bassford.  Further  information 
can  be  obtained  from  William  D.  Goodchlld,  ANA-180  (609)  641-8200,  extension 
2396. 
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TEST  NUMBER _ ^  #2  #3  #6  #7  #8  #9  #10  #12  #14  #15 

part  1  1  I  2  1  I  2  1  I  2  I  3  1  1  1  1I2  1I21I2  1I2 
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TABLE  2.  ARSR-1 

WEATHER  DATA 

Date 

Tape  Number 

Test  Number 

Time 

9ium 

1 

6 

1111 

14 

1122 

15 

1126 

16 

1130 

2 

1132 

9113111 

2 

6 

1207 

3 

1211 

2 

1217 

3 

1221 

3 

1226 

9/22/77 

1 

1 

1530 

2 

1534 

3 

1539 

1 

1544 

2 

1548 

3 

1552 

9/22/77 

2 

10 

1636 

10 

1640 

10 

1644 

10 

1648 

10 

1653 

10 

1657 

10 

1700 

9/23/77 

1 

1 

1335 

6 

1339 

12 

1344 

17 

1348 

17 

1352 

18 

1356 

^nzni 


2 


1 
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TABLE  2.  ARSR-1  WEATHER  DATA  (CONTINUED) 


Date 

Tape  Number 

Test  Number 

Time 

9/23/77 

3 

6 

1458 

1 

1502 

2 

1506 

3 

1511 

1 

1515 

6 

1518 

10 

1522 

9/23/77 

4 

6 

1546 

1 

1550 

12 

1554 

18 

1558 

1 

1602 

6 

1604 

1 

1608 

9/28/77 

1 

7 

931 

8 

932 

9 

933 

1 

937 

1.5 

937 

12 

939 

1 

941 

1 

1.5 

943 

1 

945 

12 

948 

10 

949 

1 

953 

1.5 

954 

9t2^m 


2 


1 
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TABLE  2.  ARSR-1  WEATHER  DATA  (CONTINUED) 


Date 

Tape  Number 

Test  Number 

Time 

9/28/77 

3 

Calibration 

- 

9/28/77 

4 

1 

1506 

7 

1509 

8 

1512 

9 

1515 

6 

1518 

1.5 

1523 

12 

1327 

9/28/77 

5 

1 

1554 

7 

1557 

8 

1600 

9 

1603 

1.5 

1606 

« 

12 

1610 

9/28/77 

6 

1 

1634 

7 

1637 

8 

1640 

15 

1643 

1 

1646 

1.5 

1649 

7 

1653 

1 

1656 

6 

1700 

9/28/77 

7 

1 

2100 

2 

2103 

1.5 

2106 

10 

2109 

14 

2113 

1 

2116 

1.5 

2119 

12 

2122 

9/28/77 

8 

1 

2142 

14 

2145 

2 

2147 

1.5 

2150 

10 

2152 

1 

2156 

3 

2159 

14 

2201 

1.5 

2204 

12 

2207 
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XABli  2,  ARSR-1  WEATHER  DATA  (CONTINUED) 


Tape  Number  Test  Number  Time 

9/28/77  9  7  2220 

8  2223 

9  2226 

6  2229 

15  2232 

12  2235 

1  2239 

1.5  2242 

1  2246 

9/28/77  10  1  2332 

7  2339 

8  2341 

9  2344 

1  2346 

1.5  2351 

1  2355 

10/4/77  1  1  1247 

7  1251 

8  1254 

9  1257 

1  1300 

1.5  1303 

12  1306 

12  1309 

7  1312 

-10/4/77  2  8  1320 

9  1323 

1  1326 

7  1329 

8  1331 

9  1334 

1  1337 

1.5  1340 


10/4/77 


3 


1 


^  *5  r  r 
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TABLE  2.  ARSR-1  WEATHER  DATA  (CONTINUED) 


Date 

Tape  Number 

Test  Number 

Time 

10/7/77 

1 

7 

1901 

8 

1903 

9 

1905 

1 

1907 

1.5 

1910 

7 

1912 

8 

1914 

9 

1916 

7 

1918 

1 

1920 

10/7/77 

2 

7 

1935 

8 

1937 

9 

1939 

1 

1941 

1.5 

1943 

12 

1947 

7 

1950 

8 

1952 

9 

1954 

1 

1956 

7 

1958 

10/7/77 

3 

7 

2007 

8 

2011 

9 

2014 

1 

2017 

7 

2020 

7 

2023 

8 

2026 

9 

2029 

1 

2032 

10/7/77 


1 


20A9 
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TABLE  2.  ARSR-1  WEATHER  DATA  (CONTINUED) 


Date 

Tape  Number 

Test  Number 

Time 

10/7/77 

5 

7 

2131 

S 

2135 

9 

2138 

1 

2141 

1.5 

2144 

12 

2147 

8 

2154 

12 
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TABLE  3.  DATA  PROVIDED  TO  APPLIED  PHYSICS  LABORATORY 


VQR 

ARSR-1 

WSR-57 

Test 

Window 

(Number) 

Date 

Time 

Range  Azimuth 

Range  Azimuth 

Number 

1 

9/28/77 

2142 

31  miles  266® 

40  miles  274® 

1 

2 

9/28/77 

2142 

98  miles  34® 

100  miles  28® 

1 

3 

10/7/77 

1936 

58  miles  18® 

64  miles  10® 

7 

1942 

59  miles  19® 

64  miles  10® 

1 

A 

9/28/77 

2109 

45  miles  276® 

55  miles  283® 

10 

5 

lO/A/77 

1326 

34  miles  238® 

40  miles  254® 

1 

6 

9/28/77 

2220 

6  miles  270® 

17  miles  300® 

7 

7 

10/4/77 

1329 

6  miles  270® 

17  miles  299® 

7 

13 


FIGURE  1.  ARSR-1  DATA  COLLECTION  CONFIGURATION 


.  s  s 
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78-49-3-LR 


FIGURE  3.  NAFEC  DATA  REDUCTION  CONFIGURATION 


THE  EFFECT  OF  GROUND  REFLECTIONS 
ON  WEATHER  REFLECTIVITY  MEASUREMENTS 

INTRODUCTION 

During  1911,  an  experiment  was  run  in  Oklahoma  where 
a  Federal  Aviation  Administration  (FAA)  ARSR-ID  was  used  to 
measure  storm  reflectivity  simultaneously  with  a  National  Weather 
Service  (NX'JS)  WSR-57,  as  described  previously  in  the  body  of  this 
report.  The  purpose  of  the  experiment  was  to  evaluate  the 
ARF.R-lD's  capability  to  measure  weather.  The  analysis  of  the 
experimental  data  indicated  that  the  ARSR-lD  reflectivity  values 
were,  on  the  average,  4  dBZ  higher  than  the  WSR-57  for  six 
storms  located  at  ranges  between  55  (29.7)  and  180  km  (97  nmi)  . 

In  searching  for  the  source  of  this  difference,  Ken  Coonley  of 
the  FAA  suggested  that  the  ARSR-ID  signal  may  have  been  enhanced 
due  to  reflected  energy  from  the  surface.  Reflected  energy  w  s 
a  consideration  since  the  ARSB-ID  was  positioned  with  maximum 
gain  at  3/4*  elevation  during  the  experiment.  This  positioning 
directs  a  considerable  amount  of  energy  to  the  ground. 

The  analysis  reported  herein  was  performed  to  determine 
whether  ground  reflected  power  could  account  for  a  portion  of  the 
4  dBZ  enhancement  of  the  ARSR-ID  measured  values. 

THEORY 

The  electric  field  illuminating  an  object  at  some  slant 
range  R  from  a  transmitter  can  be  expressed  as  (see  Figure  D-1) 


where 

=  transmitted  gain 
R  =  direct  path  length 
=  transmitted  power 

K  =  a  constant  which  depends  or.  the  radiating  source 
c  =  speed  of  light 
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If  the  transmitter  has  a  beam  shape  such  that  a  portion  of  the 
transmitted  beam  illuminates  the  surface,  then  the  power  incident 
at  the  target  is  the  sum  of  the  energy  transmitted  along  R  and 
that  transmitted  along  (i.e.,  reflected  from  the  ground). 

This  combined  field  is 


E 


+  4R' 


cos 


(D.  2) 


where 

=  gain  transmitted  along  indirect  path 
p  *  surface  reflection  coefficient 
D  ■=  divergence  factor  (to  be  discussed) 

R  “  direct  path  length 
Rj^  «  indirect  path  length 

e  .  aijLl  .  4,  .  .  4,  (c. 

where 

X  «  wavelength 
ij/  s  phase  angle 

L  *  path  length  difference  for  a  curved  earth 
f  »  frequency 

In  order  to  determine  A  for  a  curved  earth  we  use  the 
fol lowing  approach  from  References  D-1  and  D-2. 


Let  u  =  and  v  =  - 

*^1  /5inq; 

where 

a  =  a  modified  earth's  radius  (i.e.,  4/3  earth  was 
used  in  this  analysis) 
h,  «  height  of  the  antenna 


Page  Four-D 

=  height  of  some  point  in  a  storm 
d  =  ground  range  to  the  storm 

Consider  a  parameter  S  which  is  defined  as 

S  =  ^  (see  Figure  D-1  for  d^)  (D.4) 


This  parameter  S  is  the  root  of  the  equation 


- 1 


(D.5) 


Solving  this  equation  for  S  we  can  find  the  path 


difference 


(D.6) 


S  is  also  used  to  compute  a  divergence  factor.  When  a 
ray  illuminates  the  surface  the  beam  may  spread  depending  on  the 
geometry,  causing  a  decrease  in  the  intensity  of  energy  reflected 
from  the  surface.  We  must,  therefore,  multiply  the  reflection 
coefficient,  p ,  by  a  parameter  D.  D  is  given  by  the  equation 


Having  A  and  D,  the  field  in  Eq.  (D.2)  can  be  computed. 
The  pattern  of  energy  that  would  be  incident  on  a  storm  at  40  km 
(21.6  nni)  is  shown  in  Figure  D-2.  Plotted  in  that  figure  is 
the  quartity  under  the  square  root  in  Eq.  (D.2).  This  field  is 
called  the  interference  pattern.  One  should  observe  that  at 
certain  altitudes  the  direct  and  indirect  paths  add  to  give  twice 
the  amplitude  and  in  other  cases  they  subtract  co  create  a  null 


Range:  40  km 
Antenna  height:  10  m 
No  divergence  factor  included 
Phase  angle:  180*^ 

Reflection  coefficient:  1 
Constant  gain  assumed 


T 


0  200  400  600  800  1000 

Height  (m) 


Relative  field  strength  of  the  sum  of 
direct  and  indirect  rays  illuminating 
a  storm  at  40  km  for  the  ARSR-lD 
positioned  with  maximum  gain  at  3/4® 
elevation. 


rigure  D-2. 
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field  at  that  altitude.  Figure  D-2  is  meant  to  show  the  lobing 
pattej n  and  assumes  equal  gain  along  both  paths  with  no  divergence. 

We  wish  to  investigate  the  effect  this  interference 
pattern  may  have  on  a  storm  positioned  at  various  ranges. 

TOTAL  POWER 

Let  us  position  a  storm  at  ground  range,  d,  and  examine 
the  total  power  received  from  a  slice  of  the  storm  such  as  that 
pictui ed  in  Figure  D-3.  In  the  computation  we  divide  the  storm 
cell  into  height  intervals  Ah  and  compute  the  power  received  from 
each  small  pulse  volume.  The  power  received  from  only  the  direct 
path  ran  be  expressed  as 


K  Z(h)  G_(h)  G^(h) 

P  (h)  ^  =  -E - 1 - - 

r  direct  „2 


(D.8) 


where 


P  TT^  e  CT  |Kt^  L 

K  =  — = - - - - 

P  512 

and  P  ’..s  the  transmitted  power,  9  is  the  horizontal  beamwidth, 

^  2 

ct/2  is  the  range  resolution,  |k|  is  a  function  of  the  complex 
index  of  refraction  and  L^  is  the  radar  losses. 

The  power  received  from  the  combined  direct  and  indirect 
path  is  Eq .  (D.8)  multiplied  by  the  square  of  the  quantity  under 
the  sciuare  root  sign  in  Eq.  (D.2). 

'■  "’’total  -  S'”’ 

R 

(D.9) 
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Figure  D-3.  Sketch  showing  direct  and  indirect  rays 

illuminating  a  storm  cell  at  ground  range  d. 
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Eq .  (D.9)  represents  the  total  power  that  would  be 
received  from  the  two  paths  and  includes  the  effect  of  the 
reflected  ray.  This  effect  is  controlled  in  the  equation  by 
the  quantity  in  brackets. 

DISCUSSION 

Eqs.  (D.8)  and  (D.9)  were  computed  for  a  profile  by 
positioning  the  profile  at  various  ranges.  Since  the  interest 
in  this  analysis  lies  in  the  difference  between  the  total  and 
direct  power,  the  shape  and  magnitude  of  the  Z  profile  is  of  no 
importance  as  will  be  seen  later. 

Throughout  the  analysis  a  ground  reflection  coefficient 
p  =  1  has  been  used  along  with  a  phase  angle  of  180°.  The  p  value 
represents  a  worst  case,  which  means  that  the  difference  between 
the  total  power  and  the  direct  power  will  not  be  greater  for  any 
other  reflection  coefficient.  Since  an  integrated  power  effect 
is  being  sought  here,  the  phase  angle  simply  acts  as  a  reference 
angle  for  the  lobing  structure.  A  value  of  180°  was  used  through¬ 
out  the  analysis.  It  appears  that  at  the  ARSR-ID  frequency  and 
horizontal  polarization  a  reflection  coefficient  close  to  1  and 
phase  angle  of  180°  may  be  expected. 

FiguresC-4throughD-6show  plots  of 


S  ^r^h^  direct  ^  ^r^^^otal 

h  h 

at  three  ranges  40,  80,  and  120  km.  These  curves  were  computed 

for  \h  =  10  meters  and  to  the  altitude  for  which  the  difference 

-  53^,-  became  constant, 

total  direct 

Consider  Figure  D-4,  the  difference  in  the  two  curves 
is  small,  below  130  m  (426  ft) .  The  difference  then  increases 
up  to  210  meters  (689  ft) .  Note  that  the  summed  power  for  the 


Z  power  (dBm) 
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Figure  D-4 .  Integrated  power  at  40  kro  range  versus  Z 
height  for  ARSR-ID  radar  positioned  with 
maximum  gain  at  3/4“  elevation.  Direct 
and  total  {direct  plus  indirect)  paths 
are  shown. 


Height  (km) 


Integrated  power  at  120  km  range  versus  Z 
height  for  ARSR-lD  radar  positioned  with 
maximum  gain  at  3/4®  elevation.  Direct 
and  total  (direct  plus  indirect)  paths 
are  shown. 
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total  has  an  oscillation  in  it.  This  oscillation  is  due  to  the 
interference  pattern.  That  is,  when  the  altitude  interval  is 
in  a  null  no  power  is  added  to  the  summation.  When  a  peak  is 
within  the  altitude  interval,  power  is  added  to  the  summation. 
Therefore,  the  flattened  portions  of  the  total  power  curve 
correspond  to  the  nulls  of  the  interference  pattern. 

As  the  power  is  summed  higher  in  altitude,  the  indirect 
path  no  longer  affects  the  power  and  the  difference  between  the 
two  becomes  nearly  constant.  The  same  general  trend  is  seen  for 
80  (43)  and  120  km  (65  nmi)  except  that  the  maximum  difference  is 
greater  and  occurs  at  different  altitudes. 

The  ratio  of  the  J^total  power  to  ]^direct  power  is 


P  ' 


jd-R')  ^  +  4R'  cos^ 


This  ratio  is  independent  of  reflectivity,  Z.  It  is, 
however,  dependent  on  range  and  storm  height  since  these  parameters 
are  iricluded  in  the  computation  of  the  quantities  in  brackets. 

The  quantity  10  log  P',  power  difference,  is  shown  plotted  in 
Figure  D-7  for  four  ground  ranges.  The  curves  for  each  range 
start  at  the  radar  horizon. 

At  40  km  (21.6  nmi),  if  a  storm  is  present  which  is 
greater  than  2800  m  in  height,  the  curve  indicates  that  the  power 
received  from  the  storm  will  be  .45  dB  greater  than  it  would  be 
without  reflected  signals.  The  curve  reaches  its  maximum  at  210  m 
(689  ft)  and  decreases  with  oscillations  above  that  altitude. 

This  maximum  at  40  km  (21.6  nmi)  range  would  never  be  measured  in 
the  received  power  unless  the  storm  height  were  less  than  210  m. 

(689  ft)  which  is  highly  unlikely.  Considering  the  fact  that 
most  severe  storms  are  greater  than  5  k.  in  height,  the  curves 
indicate  that  the  potential  increase  in  received  power  would  range 
between  .5  and  2.4  dB  for  ranges  up  to  200  km  (108  nmi) . 


lOLOGP'(H)  (dB) 
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The  difference  plotted  in  Figure  D-7  increases  with 
range.  This  effect  is  due  to  the  fact  that  as  the  range  increase  ; 
the  depression  angle  of  the  reflected  ray  becomes  smaller  and  thus 
the  portion  of  the  gain  pattei  n  incident  on  the  rain  would  be 
clo£;er  to  the  maximum  gain.  The  result  of  this  is  that  the  quantity 
P'  3.n  Eq .  (D.2)  is  closer  to  1  and  the  difference  between  the 
direct  and  total  integrated  power  at  5  km  height  increases. 

Figures  D-4  through  D-7  were  computed  for  the  ARSR-ID 
antf:nna.  It  is  expected  that  other  ARSR  radars  with  narrower 
bear.widths  would  show  less  increase  in  signal,  with  the  magni¬ 
tude  depending  upon  antenna  tilt. 

For  comparison,  the  antenna  pattern  for  the  WSP-57  was 
run  through  the  program  with  its  maximum  gain  positioned  at  1/2° 
and  0°.  These  results  ere  shown  in  Figure  D-8  and  indicate  tha^- 
or  1/2°  positioning  of  the  me  ximum  there  would  be  a  small  incr(;ase 
in  power  for  profiles  at  40  (  1.6)  or  160  km  (86  nmi) .  However, 
witli  maximum  gain  at  0°,  a  3  (B  increase  in  received  power  would 
be  experienced.  This  is  probably  of  little  consequence  in  WSR 
operational  modes  since  0°  positioning  is  probably  never  used. 
However,  in  research  radar  situations,  this  increase  may  be  of 
importance  when  operatirg  near  the  horizon. 

CONCLUSIONS 

As  pointed  out  earlier,  this  analysis  was  performed  to 
determine  whether  the  4  dB  discrepancy  in  tlie  Oklahoma  experiment 
could  be  explained  by  ground  'eflections.  The  discrepancy  can 
only  be  partially  explaned  in  this  way  for  two  reasons:  (1)  The 

increase  ir  received  power  due  to  ground  reflected  rays  is  range 
dependent  and  no  range  <  ependency  v/as  evident  in  the  Oklahoma, 
data.  Il'wevcr,  it  is  pcssible  that  a  1  or  1 .  5  dB  range  dependence 
would  not  lx  rocogni  zap  e.  (2)  All  the  Oklahom.i  data  were  taken 
at  lo;  than  180  km  (87  nmj  )  ranee  and  the  largest  increase  cal¬ 
culated  (  rom  the  analysis  is  2.4  dB  at  200  km,  (10  8  r.ni)  .  It  is 
possible  that  of  the  4  dB  di'crejancy,  approximately  2  dB  could 
be  explained  for  the  ARSR  in  this  manner. 


LOG 
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Althouch  this  analysis  was  performed  for  a  specific 
reason,  it  should  be  heeded  by  those  operating  radar  with  maximum 
gain  near  the  hoiizon  when  observing  diffuse  targets.  The 
potential  should  be  considered.  In  some  cases  where  the  diffuse 
targets  are  weak  this  enhancement  may  aic  in  detection  capability 
In  cases  where  high  accuracy  is  required,  a  potential  error  may 
exist. 
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